














Comparison of Four National Radium Standards 


Part 1. Experimental Procedures and Results 


T. I. Davenport, W. B. Mann, C. C. McCraven, and C. C. Smith 


Part 2. Statistical Procedures and Survey 


W. S. Connor and W. J. Youden 


Part 1 


Che two United States primary radium standards have been compared with the British 
primary radium standard and the Canadian national radium standard (1) by an ionization 
method, using the NBS standard electroscope, 2) calorimetrically, using the Peltier-cooling 
radiation balance, (3) by Miller counter, and (4) using a scintillation 


means of a Geiger 


counter 
methods should, 
radium standard 
source self-absorption 


1. Introduction 


During January and February 1954 the British 
primary radium standard and the Canadian national 
radium standard were at the National Bureau of 
Standards for the purpose of comparing these 
standards with the two United States primary 
radium standards at the Bureau. The intercom- 
parisons were conducted over a period of 12 days 
and were made as exhaustive as possible, using the 
NBS electroscope, a Peltier radiation balance, and 
Geiger-Mialler and seintillation counters 


2. Historical Background 


In August 1911 Mme. Pierre Curie prepared, in 
Paris, a primary radium standard consisting of 
21.99 mg of the pure anhydrous radium chloride that 
had been used to determine the atomic weight of 
radium as 226.0. This 21.99 mg of radium chloride 
was sealed into a glass tube 32 mm long, having an 
internal diameter of 1.45 mm and a wall thickness 
of 0.27 mm 

At the same time Professor Otto Hénigsehmid, in 
Vienna, made three radium-standard preparations 
from very pure radium chloride consisting of 10.11, 
31.17, and 40.43 mg of radium chloride sealed in 
glass tubes about 32 mm long, having internal 
diameters of 3.0 mm and wall thicknesses of 0.27 
mm, each tube having a platinum wire sealed in one 
end. This wire was presumably to prevent the 
accumulation of static charge within the tubes. 
The purity of the radium chloride was defined by a 
radium atomic-weight determination, resulting in a 
value of 225.97. Of these the 31.17-mg preparation 
was chosen as a secondary standard. Mme. Curie’s 
21.99-mg primary standard and Professor Hénig- 
schmid’s secondary standard are generally and re- 
spectively referred to as the 1911 Paris and Vienna 
radium standards. 


317164 54 


Where there is little or no difference in gamma-ray source self-absorption, the four 
and in fact do, give good agreement 
the difference in the results obtained is an 


In the case of the Canadian national 
indication of a difference in 


In 1934, after 25 vears had elapsed, some concern 
was felt lest the Paris primary standard, together 
with a number of secondary radium. standards, 
might explode on account of the accumulation of 
helium and chlorine and possible devitrification of the 
containing tubes. Hdénigsechmid was at that time 
carrying out, in Munich, a further determination of 
the atomic weight of radium, and accordingly the 
International Radium Standards Commission asked 
him to prepare new standards, using the same salt as 
for the atomic-weight determination 

For his atomic-weight determination, which was 
carried out in the early part of 1934, Hénigsehmid 
used approximately 4 g of radium chloride, containing 
3g of radium element, that had been placed at his 
disposal by the Union Miniére du Haut Katanga 
This salt was purilied by Hénigsehmid to a poimt 
where spectroscopic analysis by Gerlach showed a 
maximum of 0.002 to 0.003 pereent of barium atoms 
\ value was obtained for the atomic weight of radium 
equal to 226.05, which is currently accepted 

Hoénigschmid then used some 817 mg of this highly 
puriiied anhydrous radium chloride to prepare 20 
new standards of radium. Exactly who asked him 
to do this is not now quite clear. According to Mile 
Chamié [1]! the International Radium Standards 
Commission, at the of Stefan Mever, 
‘entrusted Mr. O. Honigschmid with the prepara 
tion of 20 standards, using the salt he had purified 
and used in measuring the atomic weight of radium.” 
According to Hénigschmid himself, however, in a 
paper [2] presented after his death by Stefan Meyer, 
the 20 standards were prepared “at the wish of the 
Belgian radium company.” These two versions are, 
however, not irreconcilible if one assumes that the 
suggestion of the Belgian company was made known 
to the International Radium Standards Commission, 
which then gave it its official sanction, 
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The 20 new Hoénigschmid standards were sealed 
into glass tubes on -Poae 2, 1934, the glass tubing 
being similar to that used to seal the 1911 Vienna 
standard and having an internal diameter of 3.0 mm 
and «a wall thickness of 0.27 mm \ platinum wire 
wis sealed into the end of each standard 

One of the new Hoénigsehmid standards that was 
12 mm long and contained 22.23 mg of radium chlo 
ride was selected as the new international standard, 
and its value was carefully compared with the 1911 
Paris standard by gamma-ray measurements over a 
period of 4 vears | 1] The Hoénigsehmid reference 
number for this standard is 5430. Hdénigsehmid 
states [2] that the error of a single weighing was not 
more than 0.02 mg. The gamma-ray comparison 
with the 1911 Paris standard showed a discrepancy, 
however, of 0.2 percent, corresponding to a weight 
of 22.27 my as of June 2, 1934 

The first United States radium standard was 
brought to America in 1913 by Mme. Curie. This 
souree contained 20.28 mg of radium chloride and 
was designated by the International Radium Stand 
ards Commission number LV (Vienna No. 6) 

In 1936 two of the twenty Hoénigschmid prepara 
tions were acquired as the United States primary 
They are each designated by two 
numbers, namely, 5437, XIV and 5440, X\ The 
arabic numerals are those given by HlOnigsehmid, 
and the roman numerals are those assigned by the 
International Radium Standards Commission and 
imply that the standards have undergone gamma-ray 
comparison with the 1911 Paris and Vienna stand 
ards. The lengths of these two United States 
standards are 36 and 37 mm, and they contained 
50.22 and 26.86 mg, respectively, of radium chloride 
as weighed by Hénigschmid on June 2, 1984. These 
weights correspond to 38.23 and 20.45 mg of radium 
element. The weights derived from a comparison 
with the Paris and Vienna 1911 standards corre 
sponded, however, to only 38.13 and 20.38 mg, re 
spectively, of radium element, as of June 1934 

The British primary radium standard is designated 
by one number only, namely, 5432. It is solely a 
standard by weight and was not compared with the 
1911 Paris and Vienna standards. It is, however, 
one of the original Hénigsehmid preparations sealed 
on June 2, 1934. Its length is 38.8 mm, and its salt 
content corresponds to 15.60 mg of radium element, 
as of that date. This standard replaced the first 
British radium standard, which had been in the 
custody of the National Physical Laboratory since 
1913. This earlier standard was designated by the 
International Radium Standards Commission num- 
ber IIL (Vienna No. 3 

The United States and British primary radium 
standards, as can be seen from figure 1, have low 
ratios of volume of salt to volume of tube. It is 
therefore to be expected that with the standards in 
a horizontal position and the grains of radium chlo- 
ride distributed evenly along the tube their gamma- 
ray source self-absorption would be very nearly the 
same 

The Canadian national radium standard is however 
shorter and of smaller diameter than the Hénig- 


radium standards 


schmid-preparations, and it is tightly packed (fig. | 
It was sealed in June 1930 by the Union Minidre 
Haut Katanga, its contents, and that of six oth: 
sources in the custody of the National Researe! 
Council, having been taken from two tubes of radiyy 
chloride that had been prepared by the Union 
Miniére in June 1924. Its weight was derived } 
gamma-ray comparison in 1933, in Paris and Vienn»s 
with the 1911 standards, and it is designated by th 
number XIII It is understood that no corrector 
for possible differences in self-absorption were mad 
in these gamma-ray comparisons. Its length is 10 
mim, its internal diameter 1.5 mm, and its salt conten: 
corresponds, according to the gamma-ray compariso: 
with the 1911 
radium element, as of June 1934. Information o) 
all four national standards is summarized in table | 
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In view of the uncertainties that exist and the 
differences between the Hoénigschmid weights and 
their weights as derived by comparison with the 


1911 Paris and Vienna standards [3], it has recently 
been suggested that new radium standards be pre 
pared from about | g remaining of Hénigschmid’s 
original “atomic-weight”’ material. Another possi- 
bility lay in a recheck of the present standards 
With this end in view, United States primary radium 
standard 5440, XV was taken to the United Kingdom 
in the summer of 1952 and to Canada in the autumn 
of the same year. At the National Physical Labora 
torv. in Teddington, and at the National Research 
Council Laboratories, in Ottawa, it was compared, 
by gamma radiation, with the British primary 
radium standard, 5432 |4] and the Canadian national 
radium standard, NII [5]. The results obtained by 
these laboratories are discussed later in) connection 
with the data given in table 3 

The question also arises as to what is desired in a 
radium standard Im order to derive the mass of 
any radium preparation me terins of the standard by 
gamma-ray measurements it 1s know 
both the absorption of the containers of the prepara 
tion and standard and also the self-absorption of 
the radium salts themselves. In NBS certificates the 
results are stated in terms of milligrams of radium 
when contained in a Thiiringen glass tube having a 
wall thickness of 0.27 mm, together with an empirical 
absorption correction for the container in question 
Onlv calorimetric measurements can give the ratios 
of the true radium contents, 
tion but in this case it is necessary to know the date 
of sealing of the preparation in order that correction 
made for the growth of polonium. A small 
fraction of the gamma-ray energy is absorbed and 
measured by the calorimeter, but any difference in 
absorption between two sources will represent only 
a small correction to the already small contribution 


(about 7% to the 


hecessary to 


irrespective of absorp 


may be 
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total energy emission, 
3. Measurements With the NBS Standard 
Electroscope 
The NBS standard clectroscope 16] and Measuring 
were without modification, for this 


comparison of four national radium standards. The 
ionization chamber consists of a 10-em cube free-air 
| 


svstom used, 


volume, with walls made of em of lead and a 
‘cm aluminum inner lining. A gold leaf is) sus 
pended near the center of the chamber, A> 10-u 


quartz fiber at the free end of the leaf provides a 
fine line for projection. The fiber image is magnified 
approximately 100) times and projected onto a 
metric seale. The discharge of the electroscope is 
measured by timing the transit of the image between 
two fixed points on the scale 6 em apart 

The source indexing system consists of a V-shaped 
trough of “%e-in. Lucite on an aluminum stand 
The stand can be moved along a line perpendicular 
to the face of the ionization chamber or rotated 
about its own vertical axis Preparations are 
centered in the trough opposite the center of the 


chamber, so that measurements are made perpen- 
dieular to the axes of symmetry of the preparations 

The four standards were measured relative to each 
other by comparison of each of the six possible com- 
binations of pairs. Independent measurements were 
made on each pair by each of three different observers 
at source distances of 66.5 em and 74.1 em from the 
chamber. The entire series of measurements was 
repeated twice 

The following procedure was adopted for compar 
ing each pair of standards 

1. The trough was placed at the distance selected 
and parallel to the chamber face 

» A standard was held horizontally and tapped 
lightly until the salt was distributed uniformly along 
the length of the capsule, as in figure 2 

}. The standard was placed in the trough and cen 
tered 

+. Three observations of the discharge time were 
made and recorded 


5. The trough was rotated 180 degrees, and three 
more observations were made 

6. Procedures | to 5 were repeated with the second 
standard of the pan 

7. Procedures 1 to 6 were repeated for both mem 


bers of the pair at the second distance from the 
clectroscope 
s 
| 
cr 
- 
at 
A 
ricvure 2 Four national radium standards, with the qrain 
of salt in the three Hénigschmid standards distributed a ong 


the length of the tube 
1, American: A, British: C. Canadian: and DP, Ame 


4. Comparison by Geiger-Muller Counter 


The Geiger Miller counter used for this compari 
son was a neon-halogen-filled tube. The tube itself 
was surrounded by a sheath of lead % inch thick so 
that the soft gamma rays, the spectrum of which 
might be varied by source absorption to a greater 
extent than that of the higher-energy gamma rays, 
would not be counted. The resolving time of the 
counter was determined by the two-source method 


to be 211 psec 5 percent The correction for re 
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solving time applied to the data ranged from 1.1 to 
2.7 percent 

The source holder of the NBS standard electro 
scope was used to position each standard in turn in 
these mensurements, and the standards were tapped 
so that, in the ease of the more loosely packed 
Hénigschmid standards, the grains would be dis 
tributed uniformly along the tube 

In order to elimimate any possible effects due to 
drift a series of measurements was carried out on each 
pair of international standards. Thus, in the com 
parison of A and B, measurements were carried out 
with A and B arranged in “packages” in the following 
order: A, B, B, A; B, A, A, B; A, B, B, A; and 
finally, B, A, A, Bo Similar package measurements 
were made on each of the other five pairings of the 
four international standards 

\ total of about 80.000 counts was taken on each 
of the 16 members of the 4 packages comprising a 
pair comparison. Thus in the comparison of A and 
Bea total of some 640,000 counts were made with A 
in position and 640,000 with B 


5. Comparison by Scintillation Counter 

The seintillation counter consisted of a thallium- 
activated sodium-iodide erystal mounted on the face 
of a photomultiplier tube. The resolving time of the 
counter and amplifier was 5 use 10 percent, and 
the corrections applied to the data varied from 0.5 
to O.8 percent. The discrimimator was set to accept 
pulses corresponding to gamma-ray energies greater 
than | Mev. Thus, as for the NBS standard elec- 
troscope and the Geiger-Miller counter, the effect of 
source self-absorption of the lower-energy gamma 
ravs should not be apparent. The sodium-iodide 
crystal and photomultiplier were mounted adjacent 
to the Geiger Miller counter so that counts on each 
source could proceed concurrently with both counting 
systems. Exactly the same pairing and packaging 
order of sources as was used for the Geiger-Miller 
counter comparison was ipso facto, also used in the 
scintillation-counter measurements. The counts for 
each source in position were of the order of 400,000 
compared with 80,000 in the case of the Geiger- 
Miiller counter. 

6. Measurements With the Radiation 
Balance 


A modification of the radio-balance originally 
designed by Callendar [7] for the measurement pri- 
marily of radiant energy has recently been ieetbed 
[8], which is suitable for the measurement of the 
energy emission from radioactive materials, This 
modification of the radio-balance has been renamed 
the radiation balance, its most important feature 
being the ability to balance the energy emission from 
a radioactive source either against Peltier cooling or 
the energy emission from another radioactive source, 
or both 

None of the radiation balances constructed pre- 
viously was large enough to accommodate the large 
Hénigschmid standards, and, accordingly, a new one 
was constructed for this purpose. This balance is 


described in detail separately in this issue [9]. It 
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differed essentially from the first one, howeve: 
that its larger cups were made from gold iInstenc 
copper 


7. Radiations Measured 


The radiation from radium in equilibrium with 
its products consists of five energetic alpha partir 
groups, including that of polommum: three ny» 
groups of beta particles, the most energetic be) 
that from the transition of radium E to radium 
with a maximum energy of 1.17 Mev: and a co: 
plexity of gamma rays, the most energetic bei 
from the excited levels of radium C’ 

Three of the methods deseribed here and sed 
compare the a contents of the four nation 
radium standards were essentially gamma-ray con 
parisons With the thicknesses of lead used. o1 the 
discriminator setting, the chief contribution to thy 
Lrentritiia-ray effect would be from the ehneryvety 
radium C’ gamma ravs (above 0.6 Mev in the case of 
the clectroscope and CGemer Miller counter ana 
above | Mev in the case of the scintillation counte: 

In contrast, the radiation balance measures pri 
marily the energy emitted in) corpuscular form 
Some 93 percent of the energy produced by radium 
and its daughters down to radium D is associated 
with particulate emission, the remaining 7 percent 
of the energy produced being associated with th 
gamma radiation. The wall thickness of the gold 
cups was such as to absorb completely the most 
energetic beta particles from radium FE. Some 12 
pereent of the energy associated with the gamma 
rays is also absorbed. Of the 7 percent of the total 
energy produced that is associated with the gamma 
ray emission, another | percent (for the Canadian 
standard) or 1.5 percent (for the Hénigsehmid stand 
ards), corresponding, respectively, to 0.07 and 0.1 
pereent of the total energy produced, will be absorbed 
in the sources themselves. The difference of 0.5 
percent between the source self-absorption of the 
Canadian and Hénigsehmid standards corresponds 
therefore to a difference of only 0.035 percent of the 
total energy produced, which is negligible. Any 
smaller differences in gamma-ray source self-ab 
sorption of the three Hénigsehmid standards are also 
therefore neghgible so far as the measurements in 
the radiation balance are concerned. The alpha 
particle and beta-particle absorption is complet 
a correction must be made, however, for the growth 
of radium FE and polonium, which will not be tn 
equilibrium with the radium 


8. Results 


The results obtained with the radiation balance: 
measuring the sources singly and in every combina 
tion of pairs, are summarized in table 2, In this 
table the order of measurement is represented b 
reading from left to right and down the table 

From the results in table 2 the following best est: 
mates for the energy absorbed (in microwatts) from 
sources A, B,C, and D have been deduced: 

A B c D 
6203. 4 2569. 8 1131.0 3360.7 
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Tanure 2 Energy absorption, in microwatt 
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In table 3 are shown the complete results for the 
using the NBS. stiundard 
CGiewer-Miller countes 
counter, and radiition balance 


iX pairs of standards 
electroscope, semitillation 
In the last line of 
thy table are shown the wewht ratios for the same 
The weight of the Canadian standard 
(') is, however, only a derived weight, and for this 


IX pairs 


reason, anv ratio involving this derived weight is 
The ratios ABO AD 


the ratios of Pl dniigesc hiul's 


hown im quotation marks 
and 2) are, however 
own weighings 
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Kor comparison with these values the ratios ob 
tnined by Perry [4], using the NPL standard toniza 
tion chamber with voli leaf ele troscope, hhh toniza 
tion chamber with a Lindemann electrometer, and a 
Geiger-Miller counter for B/D) were, respectively, 
0.7609, 0.7657, and 0.7669 The result obtained for 
the gamma-ray ratio (77) by Michel [5], using the 
NRO preeision ton chamber and Lindemann elee 
trometer was 1.192 Michel, from geometrical 
considerations, then caleulated the source absorp 
tion of each standard and corrected the gamma-ray 
ratio to give a weight, or content, ratio of ¢ Dequal to 
1.185. The direct gamma-ray ratios obtained both 
by Perry and Michel are in excellent agreement with 
the results shown in table 3 

\ check on the internal consistency of the results 
shown in table 3 can be provided by assuming that 
A, Bo and Dare so much alike that there are negli 
gible differences in source absorption for high-energy 
yamma rays, and none at all in the case of the calo 
rimeter, where 93 percent of the energy absorbed is 
particulate, so that any change due to absorption 
of the 7 percent of gamma rays and secondary 
electrons would be even more negligible. A check 
can then be run on the results for A, B, and D by 
dividing the quantity characteristic of each in each 
determination by the Hénigschmid weight of each 
standard, This characteristic quantity is seale divi 
sions per second for the NBS standard ele ‘troscope, 
counts per second for the counters, and microwatts 


for the radiation balance In each ease the ehat 
acteristic quantity is as of February 1954, and the 
mass of radium element is as of June tos It is not 
necessary for this check to correct for the JO-veutr 
decay of radium as this is the sameceonstant for each 
standard 

The results of this internal precision check are 
shown in table 4, in which the figures quoted are the 
characteristic quantity, divisions, or counts pet 
second or microwatts, divided by the mass of radium 
clement present and normalized to make the “best 
nverage equal to 100.00 in each cause 
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should, im turn 
enable one to form an estimate of the preemion of 


The values of the best average 


Honiwschmid’s werht determinations, im whieh 
according to THoniwschmid himself [2], the error of 
ao single weiwhing was not more than 0.02 my \ 


statistical survey of the results was carried out with 
the cooperation of WS. Connor and W. A. Youden 
and resulted im the best estimates of the mass of 
radium element in al, Bo and J? given in table 5 
The methods adopted to arrive at these best esti 
mates, and also the best estimates given mn table 3 
for the ratios of pairs of standards, are described by 


Connor and Youden in part 2 of this paper 


Panis ’ Heat catimate n milligram of the manne of the 
lHlonig chmid radium standard as of June % 1934 
Standard | It iD) 

Honipeehmid nin ts US 1s. oO 20. AN 
Mass derived from NUS 

standard clectro cop 48, 227 15. ol 20. 446 
Ma derived from Creiget 

Muller counter tS. 2h 15. SOS mM 4d 
Mass derived from scintil 

lation counter 38. 245 15. 505 20.4404 
Ma derived from the 

radiation balances BS. 2B5 15. 608 2 ANS 


9. Mass of Radium Element in the Canadian 
National Standard 


bby comparing the calorimetric ratios yiven on 
table 3 with the “weight” ratios, itis clear that the 
derived weight of the Canadian national radium 
standard (C') is low by about 3 percent llowever, 
this does not allow for the difference in sealing date, 
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which involves a compensating polonium-growth 
correction of about 1.8 percent. By comparison of 
the “weight” ratios with the NBS standard electro 
scope ratios, it is also confirmed that no source 
self-absorption correction could have been made in 
deriving the certified weight of radium in_ the 
Canadian standard. However, from the data avail 
able it is possible to derive a value for this mass of 
radium 

The experimentally determined ratios of the energy 
absorbed in the radiation-balance cups) per unit 
mass of radium element for A, B, and /? are 164.62, 
164.73, and 164.54 ewomg, respectively. Taking 
the best average value of 164.58 sw/me of radium 
element, the mass of the radium in the Canadian 
national radium standard is found to be equal to 
25.10 my. as of June 1934, uncorrected for the growth 
of polonium or of radium E 

Using the Curie-Yovanovitch equation, as cor- 
rected for new values of the decay constants by 
Jordan |8, 10], the energy increments due to growth 
of polonium-210 in A, B&B, and DP), on the one hand, 
and in ©, on the other, are found to be equal, re 
spectively, to 12.2 and 16.2 cal g' hr“! inclusive, of 
nuclear recoil energy, the separation and sealing 
dates being, respectively, May 25, 1954, and June 
2 1934, for the Hénigsehmid standards, and June 
1924 and June 1930 for the Canadian national 
standard. The growth of radium E will contribute, 
in proportion, another 0.8 and 1.0 cal g"' hr’. Sub- 
tracting the contributions of polonium-210 and 
nuclear recoils and of radium E from the energy 
absorbed from A, B, and 7) in the radiation-balance 
cups gives a total energy absorption for all three 
sources equal to 11103.0 instead of 12223.9 ww (as 
of February 1954) 

In the case of the Canadian national standard, an 
energy production of 17.2 eal ¢*' hr-' by polonium- 
210 and radium E corresponds to 20.0 uw/mg of 
radium element, which. by a second approximation, 
is found to be equivalent to 489.7 ww/24.48, mg of 
radium element (the mass of radium as of June 
1924). The corrected energy absorption from the 
Canadian national radium standard is therefore 
3641.3 instead of 4131.0 ww, as of February 1954 
The radium content of the Canadian national 
radium standard, as of June 1934, is then obtained 
by multiplying the total weight of the Hénigsehmid 
standards (76.28 .ag as of June 1934) by the ratio 
of the corrected energy absorptions of February 
1954. This gives the result that there were 24.36 
mg of radium element in the Canadian national 
standard, as of June 1934. This value will, if any- 
thing, be on the low side, however, as some radium 
D on the walls of the original two tubes may have 
been lost on transfer when the Canadian standard 
was resealed in June 1930. In this event, the 
polonium-210 correction will have been too great 


10. Summary of Results 
As a result of this intercomparison of national 
radium standards, the ratios of the weights ascribed 
to three of them by Hénigschmid have been con- 
firmed. 





from the comparison of the two United States st; 
ards with the 1911 Paris and Vienna standards 
therefore, too low; unless it were assumed that 9 
Hénigsehmid’s mass determinations were low in 


same ratio However, this is to be discounted 
cause the Berlin standard was, by comparison \ 
the 1911 standards, found to have a greater wei 
than that determined by Hénigsehmid [1] 

Relative to the Hénigschmid weights, the Cy 
dian national radium standard is found to hay 
mass of radium element equal to 24.36 mg, wh 
indicates that no correction for difference in sou 
self-absorption was made in its comparison with | 
1911 Paris and Vienna standards. The differen 
between this value and that obtained by compariso 
with the 1911 Paris and Vienna standards (24.2% 
mg as of June 1934) would indicate a self-absorption 
correction of 0.53 pereent The absorption corre: 
tion determined by Michel [4] was 0.94 percent: thy 
difference between these two values could be 4 
measure of the loss of radium D and polonium-210 
in the transfer of June 1930 
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Part 2. Statistical Procedures and Survey 


W.S. Connor and W. J. Youden 


Thu tatistical analysis of the observations on the four national radium standards is 


discussed 


Che readings made with the elect ros ope, Geiger- Miller counter, and scintillation 


counter were adjusted by one formula, and the readings made with the radiation balance by 


a different formula 
criterion Finally, the adjusted values wer 
relation hip between the masses 


1. Introduction 


Four national radium standards were recently com- 
pared at the National Bureau of Standards, as 
deseribed in part | of this paper. The unusual 
opportunity with the presence of four 
standards in one laboratory directed attention 
to certain statistical aspects of the intercompari 
son. The experimental procedures and results are 
deseribed in part | 


nssociated 


Part 2 discusses the statistical 
analysis 

When two standards are compared, careful meas 
urements provide an estimate for the value of one 
standard in terms of the other A standard error 
mav be calculated for this estimate. A third standard 
makes possible the additional experimental evalua 
tion of each of the first two standards in terms of the 
third 

Suppose that three standards A, B, and ( are 
available. The experimental ratios a/b, b/c, e/a may 
each be determined by using exactly the procedure 
that would have been employed if just two standards 
had been available None of the measurements made 
on A in estimating a/) are used in the estimation of 
cla Additional data for A are taken to determine 
ca. There isa considerable advantage in this method 
because the precision of the comparison is improved 
by alternating the readings on the two standards 
under comparison. This alternation reduces the 
effects of drift in the instruments and changes in the 
environment. As soon as the ratios a/b, b/e, c/a have 
been determined there is a simple test for the con- 
sistency of the three ratios. The product of the three 
ratios should be unity. The diserepaney between 
this product and unity provides a measure of the 
errors in these ratios 

\ similar consistency criterion was applied to the 
six ratios determined by the electroscope, Geiger- 
Miller counter, and scintillation counter. Because a 
different statistical treatment was required for the 
measurements made with the radiation balance, 
those measurements are discussed separately 

The last section describes how the masses of the 
standards were used further to improve the estimates 
of the standards 


2. Comparison of the Standards by Means 
of Electroscope, Geiger-Miiller Counter, 
and Scintillation Counter 


Using these methods, environmental conditions 
common to paired measurements introduce a common 
multiplicative error in the measurements 


It is ad- 


In each case the adjusted values 
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of the standards satisfy a consistency 


improved by making use of the proportional 
and the radioactive effects of the standards 


Vantageous to express the results of paired measure 
ments as ratios to eliminate this error 

There were four standards, A, B, Cy and). There 
fore, the following SIX ratios could be determined 
experimentally 


a h ae ad d h ( 4 d ‘ d 


These provide opportunities to test the consistency 
of the data Kor example, the products 


ab «bie e/a 
ahxbdxdia 
aexecdx<dia 
hexedxdb 


should all be equal to unity. The discrepancies 
between these products and 1.0000 reveal the errors 
of the measurements. It is proper to make use of 
the information that the products should be exactly 
equal to one. The measured ratios may be adjusted 
by a least-squares technique to obtain new ratios 


a a “ “ ‘ 
A/B, A/C, ete., which do in facet multiply out to 
unity for all combinations that should give unity. 


This includes not only three factor combinations 
such as 

o “ “ 7. A . 

A/BX BICXCI/A 


but also four factor products 


Nhs Bibs Oops DA. 


The adjustment formula used on the data shown 
in table Lis of the form 


V5) (<5) 


where the lower case letters indicate the measured 
ratios.’ The adjusted values (see table 3 in part | 


——_—__—_——_. 


SUP | > 


lhisadjustment formula is related to the adjustment formula for the difference 
between the estimates of two treatment effects in a balanced incomplete block 
(BIB) design, see K. L. Anderson and T. A. Bancroft, Statistical theory in re 
search, p. 252 (MeGraw-Hill Book Co., Inc., New York, N. Y., 1982 Since the 
two measurements in a pair, as a and or cand d, are subject to a common multi 
plicative error, the logarithms of the two measurements in the pair are subject to 
4s common additive error, Hence, the BIB design formula applies for the differ 
ence between the logarithms of the adjusted values, as bg A~loe HF, and by 
taking antilogarithms, the above formula ts obtained 


have the property that 


The observed values do not meet this consistency 
requirement. The reconciliation among the results 
effected by the above least-squares technique intro- 
duces each standard symmetrically in the computa- 
tion pattern and does not single out any one standard 
as a superstandard After the relative values have 
been established standard may be given 
value, whereupon all other standards 
determined without changing the relative values 


one an 


avreed ure 


Panis | / rperrmenta f ; / ifios of fou landard 
Ni 1 
hl t 1™ tis ‘ s teh wd 101s 
Cheer MOller ¢ ninter >ayu j s 1 sys4 4450 “Tso 1 10 
tilhath tsi ! st) ; “ine 1 
The above least-squares adjustment has long been 


used for other comparisons. Recently, it has been 
found that certain subsets of pairs selected from all 
possible pairs lead to convenient least-squares esti 
mates” Given that a reasonably small number of 
pairs will suffice to interrelate all the standards, 
there would appear to be some chance of success for 
an international program of comparison 
properly selected subset of pairings was obtained, 
the various national standards could be tied together 
with values that would give consistent comparisons 
among the standards 


Once a 


3. Radiation-Balance Measurements 


The radiation balance used in this work was suit 
able for measuring either a proportion of the energy 
emitted by one standard or the same proportion of 
the difference in energies emitted from two standards 
This difference is determined by measurement 
The schedule of measurements included separate 
measurement on the four standards as well as the 
six possible differences between them. The pre 
cision of measurement of a difference was the same 
as the precision of measurement of a single standard 


Typical formulas for the least-squares estimates 
for the 10 quantities follow 





a 
A a lia h) b} | a ( ‘ | a d)\4 | 
A 
(L1— 8) (a—b)+1lla—e c—b)| 
Hia—d d—b)|+ha—2bb 
WJ. Youden and W. 8. Connor, Making one measurement do the work of 
two, Chem, Eng. Progr. 48, 449 (1053); and W. J. Youden and W. 8. Connor 
New experimental designs for paired observations, J. Research NBS 3&8 


19M) RP2S2 
For a discussion of the method «f least squares, see KR. L. Andersen and T. A 
Rancroft, Statistical theory in research, p. 1 Metiraw- Hill Book Co., Ine 


New York, N.Y wiv 
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The 


a a“ 
quantities. The value for (A— 8B) given by the abo 
formula will agree exactly with the result: obtain 


quantities a, 6, (a—b are measur 


7 
by subtracting the adjusted estimate B from t} 


adjusted estimate A. ‘This was not true for the 

corded values. The total amount of energy measur 
for the standards ts left unaltered by the adjustmen 
Slight shifts take place in a, 6, ¢, d, (a—b), ete, v 
achieve the The di 
crepancies between the measured quantities and th 
corresponding adjusted values afford a measure of thy 
precision of the measurements. The calculation 

It should be noted that no quan 
twice 


consistency among results 


shown in table 2 
tity 
cealed 
that ts, 
ten 


was measured The replication is con 
There are, of course, only four standard 
four quantities to be determined from th 


a 
lhis 


six degrees of freedom, available for estimating thy 


observations leaves SIX contrasts, 1. « 


standard deviation 


Papi « landard deviation n" 


(Calculation of microt 
for radiation balance 
Standard Observed Adju ted Differ Differ 
ence ene 
uu ue uu uu 
| H2S85. 2 6203. 4 8%. 2 67. 24 
R 2571.0 2560.8 1.2 1.44 
( 4127. 1 1131.0 1. 0 1h. 21 
D 3371. 6 3360. 7 io 9 118 81 
1... B 3727. | 3723. 6 3.5 12 25 
i—¢ 2164.2 2162. 4 1.8 3. 24 
i DvD 2035. 6 2032. 7 2.9 & 41 
( Kk 1561, 2 1561. 2 00 0. 00 
C—D 776. 0 770. 3 5.7 $2, 40 
DR 7RS. 6 700. 4 2.3 5. 20 
264.38 
6 Ou 


Standard deviation \ 


4. Masses of the Radium Standards 


Standards A, B, and )) were made from the samy 
supply of radium salt. The weighings were made in 
the same day by Hénigsehmid and are considered to 
have a maximum error of 0.02 mg The various 
properties of the three Hdénigschmid standards 
measured by the several methods used in this inter 
comparison are believed to be directly proportional 
All the methods 
give relative values for the standards. In addition 
the radiation balance measures the difference between 
two standards directl, Standard D) was 
arbitrarily given the value of unity and the values for 
A and B expressed relative to it. Table 3 contaim 
some of the adjusted ratios from table 3 of part 
including the ratios derived from Hdénigsehmid 
weighings 


to the masses of the standards 


any 




















| ARLI ) Value of standard, when D equa 1 tHn) 
Method 1 BR D 
Scintillation counter 1. SSO 0. 7624 1. 000 
Geiger-Miuller counter 1. SS5 7659 1. OOO 
NBS standard = elec 
troscopx 1. S70 7THG1 1. OOO 
Radiation balance 1. 873 T6A7 1. 000 
Weighing 1. SOU 7628 1. OOO 


For each method of measurement a plot may be 
made of the values of the standards against the 
corresponding masses. The resulting points should lie 
along a straight line that passes through the origin 
Let m denote the mass and r the radioactive effect 


viven by any one method of measurement in the 
relative units of table 3. If the error in r is & times 
as large as the error in m (as measured by the 


standard deviations), the slope 6 may be computed 


from the quadratic 


and r, () = «4, B,D) are the masses and the 
corresponding radioactive effects for standards A, B, 
and J) from table 3 This method of determining 
bh has the property that the sum of the squares of 


where ” 


the perpendicular distances of the points (Am,, 
from the line r= b(4m) is minimized 
‘ For a discussion of this method, see W) Edwards Demi tatistionl adjust 
t lata, Fe xe eG. 184 John Wiley & Sons, lr New York, N. ¥ NX 


mi 








In general, any particular plotted point will not be 
located exactly on the fitted line. The plotted 
points are subject to errors of observation. The 
“best” estimates of the coordinates for the point are 
taken to be the coordinates of the point on the line 
plotted pol coordinates 


nearest to the These 


; ; 
m’ and r’, are 


m and / b 


This procedure for fitting lines was followed for 
each of the lines relating the measured radioactive 
property to the mass. For the electroscope, Geiger- 
Miller counter, scintillation counter, and radiation 
balance, & was taken as 2, 4, 5, and 1, respectively 
These values correspond to the given by 
table 4 of part 1, except for the scintillation counter, 
for which 5 was used mstead of 6 


errors 


To obtain estimates for the masses, each value of 
m, Was multiplied by Hdénigsehmid’s value for /)/, 
i. «., 2045 me. These estimates are recorded in 
table 5 of part | In every case the result agreed 
with the assigned mass within the claimed weighing 
error, It is particularly interesting to observe that 
the estimates obtained from the line, using the 
radiation balance results, confirm the assigned masses 
For this line the errors in m and r were taken to be 
the same, and therefore any displacement of the 
point to bring it on the line required equal changes 
in the experimental values for mass and energy 
Mav 27, 
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A Radiation Balance for the Microcalorimetric Compari- 
son of Four National Radium Standards 


W. B. Mann 


of a radiation balance, 
three 


The design 


to accommodate Honigschmid 


element in terms of all three Honigsehmid 


u hie 


1. Introduction 


\ series of measurements has been completed in 
which the Canadian national radium standard and 
the British and United States primary radium stand 
ards were compared with each other at the National 
Physical Laboratory, Teddington, England l1],' at 
the National Research Council's Laboratory, Ottawa 
Canada [2], and at the National Bureau of Standards 
Washington, D.C. [3] 

The measurements at the National Bureau of 
Standards were made by means of the NBS standard 
1| the 


\Mliiller and se intillation counters 


clectroscope radiation balance, and Geiger 
The radiation balance already described [5] was fat 

mall to the large Hoénteschmid 
standards, which are approximately 38 mm long with 
an external diameter of mm. The Cana 
dian national radium standard is of considerably 
smaller dimensions, 10.5 mm long and 2.0-mm ex 
ternal diameter, and could be inserted into the cups 
of the first The radiation balance 
structed to accommodate these Hénigschmid stand 
ards is described here 


accommodate 


low 


about 3.5 


2. Radiation Balance 


The radiation balance is a twin microcalorimeter 
that utilizes the Peltier effect to equalize the tem 
peratures of the two cups. The utilization of Peltier 
cooling in microcalorimetry, in addition to its use by 
Callendar [6] and Hoare {7, 8], has also been sug 
gested and used by Duane |9, 10], Tian [11, 12] and 
and Calvet [13], but none of these authors has used 
the transfer method of Callendar, which eliminates 
the necessity for correction for the Joule heating in 
the Peltier couples and leads. Calvet uses differ- 
ential systems of two and four cups for the micro- 
calorimetry of slow phenomena, which might not be 
susceptible to measurement by the transfer method, 
and Swietoslawski states [14] that Duane considered 
using a compensating system of two microcalorim- 
eters, but that none of his results is to be found in 
the literature. It would appear therefore that, apart 


Figures in brackets indicate the literature references at 


he end of this paper 


a twin tmicrocalorimeter 
radium 
production of radium and its daughter products down to radium D, for unit mas 
tandards, wa 
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Peltier effect 
rate of 


utilizing the 


standards 1 described The energy 
of radium 
found to be 138.6 cal 


equal to 


from the work of Calvet” which uses both the Peltier 
and Joule effects to compensate the heats of reaction 
this method of calorimetry has not hitherto made 
any great appeal 

The cooling in the radiation balance is applied es 
sentially at a point, It is important therefore not 
to increase the of the contamers, o1 
to such an extent that temperature gradients 


Size source 
“cups” 
along them shall be so significant as to cause appre 
ciable differences in heat transfer between each of 
the cups and their surroundings. It is important 
that the thermal balance shall be effected chiefly in 
the material of the themselves In the first 
balance | the cups were coppel tubes « losed at the 
lower ends, 'j, inch long and having unternal and 
external diameters of 0.099 and O.130 inch, respec 
tively 
the accurate comparison of radium preparations con 
taining of the order of 5 or 10 meg of radium element 
By using gold cups of greater wall thickness it was 
found possible to improve the thermal mixing in the 


cups 


This balance had adequate sensitivity for 


cups and yet, at the same time, to keep the sen 
sitivity at about the same value 

A new balance was therefore constructed, and ts 
illustrated in figure L, in which the cups consisted of 
gold tubes, again closed at the lower ends, each 


18% inch in length and having internal and external 
diameters of 0.155 and 0.238 inch, respectively. The 
wall thickness of the gold cups is thus 0.083 inch as 
compared with 0.031 inch of the copper cups of the 
original balance. The same type of Peltier couples 
and thermopile junctions of chromel and constantan 
were used as before, and, with the same galvanometer 
and seale position, the sensitivity obtained was 0.21 
mm/pw at 25° C as compared with a figure of 0.35 
mm/uw at 25° C for the old balance 

Apart from enlarging the copper block containing 
the radiation-balance cups, the general details of 
construction were, as from figure 1, 
essentially the same as before, and the same auxiliary 
electrical equipment was also used 

After insertion of the standards to be measured, o1 
of the dummies of the standards, into the cups, 
silicone oil was used to fill up the cups, which were 
closed by means of stainless-steel or phosphor-bronze 
ball bearings “{¢ of an inch in diameter. 


can be seen 


An extensive bibliography ts given in [15 





rieure | Gold-cup radiation balance 


3. Theory 


If a source A dissipating energy at the rate of VW 
watts in cup | of the radiation balance be exchanged 
with a source B emitting energy at the rate of W, 
watts in cup 2, and if, at the same time, the current, 
(through the Peltier junctions, which are connected 
in series, is reversed, then 


Wy 2) d,—d,)/s, | 


where 7’ is the Peltier coefficient, d, is the galvanom- 
eter-seale deflection with A in cup 1 and B in cup 2, 
and d, is the galvanometer-scale deflection with the 
sources exchanged and the current reversed; s is the 
galvanometer — sensitivity in =o millimeters — per 
microwatlt 

The Peltier coefficient, 7’, is equal to TdT, 
where 7 is the absolute temperature, and d//dT is 
the thermoelectric power of the Peltier couple in 
microvolts per degree, In general, d/-/dT is deter- 
mined at some given temperature and must then be 
corrected slightly for fluctuations in the ambient 
temperature. In practice, d/-/d7T was expressed in 
microvolts per degree at 25° C and was then corrected 
by adding or subtracting 0.08 wv /deg C for a chromel- 
constantan couple for each degree above or below 
25° C, the value of d//dT at 25° C being of the order 
of 59 uv/deg C It is also necessary to correct s by 
the same temperature coefficient in order to allow for 
the slight change in response at different temperatures 
of the chromel-constantan couple 

The ratio of the energy dissipations of two sources 
may, however, actually be determined without 
measuring d//dT as eq (1) may be rewritten 


WW, (270+"5" or 


where «’ is now the scale sensitivity in millimeters per 
ampere per degree Kelvin, and the source in the 
second cup is a dummy, i. e., W, is equal to zero. 


As di-dT at 25° C is a constant of the apparat 
the ratio of the rates of energy production in vari 
radioactive sources can be obtained simply in te: 
of measurements of current, temperature, 
galvanometer-scale deflection. Such — ratios 
independent of the absolute value of dE/dT. 4 
therefore the uncertainty in the measurement 
dE/dT does not enter into such a determination 
relative values. Allowance must be made, howe, 
in determining such ratios for the known temperat 
coefficient of d/i-/dT in cases where Wy WwW: is det 
mined at one temperature and WY—W at anotly 

As the Peltier effects in each cup will be modifi: 
by small losses of heat to or gains from the surroun 
ings, the value of d///d7T that must be used in eq 
is not, in general, equal to the thermoelectric power of 
the Peltier couple as determined by the measurement 
of the variation of electromotive force / as a function 
of temperature fl It is necessary therefore to 
determine the effective value of dE /dT by means of « 
heating coil, together with a compensating coil 
which can be placed inside the two cups. In this 
calibration it should not be necessary exactly to 
reproduce in the heating coil the dimensions of the 
radioactive source, provided the cups are of such 
dimensions, and of material of sufficiently good 
thermal conductivity, that large differences of tem- 
perature cannot be maintained and = that cood 
thermal ‘“mixing’’ results 

In actual practice, the value of dE /dT for the 
chromel and constantan wire used in both the old 
balance [5] and that described here was found to b 
approximately 61 gv /deg C, by direct measurement 
whereas the calibrated values for both balances gave 
a value of dE /dT of about 59 uv /deg C (59.11 for the 
old [5] and 58.78 for the new, which has larger cups 


4. Calibration of the Radiation Balance 


In view of the large differences in dimensions of thy 
Hénigsehmid radium standards, on the one hand, and 
the Canadian national standard, on the other (see 
fig. 2), it was felt desirable to use two sets of calibra 
tion coils, one set simulating the Hdénigsehmid 
standards and the other simulating the smalle: 
Canadian standard, These are also shown in figur 
2. The resistances of these coils were measured 
after the calibration experiments in the radiation 
balance had been carried out in order not to impat 
the insulation of the leads to the coils, which had to 
be bared in order to attach potential leads of 38-gage 
constantan wire. On both sets of coils these poten 
tial leads were attached with an accuracy of about 

-0.1 mm to points on the coil leads approximate!) 
5 em and then 4 em from the ends of the coils that 
rested on the bottoms of the radiation-balance cups 
when the coils were inserted into them. If 2 be th 
resistance of the main coil of either set, r the resist- 
ance of the compensating coil, 7 the current im the 
coils in series, and C the current in the Peltier junc- 
tions in series, then on interchanging the coils 
between cup 1 and cup 2 and reversing the Peltier 
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ry thE 2 I ) nation l wi ; standa 
adi 437, XIN 
rT , 40, XA 
current, C, we have 
rPR—) iN NW 
' d d dk 
(27C4 = ar 
whenes 
dk L R / . 
dT 27TC+-(d,—d,)/s’ is 
As (d,—d is usually small compared with 27C 


it is sufficiently accurate, in deriving di-/d7, to use 
s as a first approximation, and to derive a slightly 
more accurate value as a second approximation 

The values obtained for the resistances of the 
calibrating coils, when compared by means of a 
potentiometer with a 10-ohm = standard coil, are 
shown in table 1. The results obtained by balancing 
the Peltier cooling against the Joule heating of the 
calibrating coils are given im table 2 

The average values for d//d7. using the results 
shown in tables 1 and 2, are 58.91 pwv/deg C for the 
“Honigschmid” coils and 58.64 ywv/deg C for the 
Canadian” coils. The average of these is 58.78 
uv deg } 

The difference between the results for the two 
coils is not considered to be outside the experimental 
error for this particular calibration. The nearer to 


‘ 


a point source of radiation to which one approxi- 
mates, the more effective will the Peltier cooling 
become, so that one might expect the Canadian coils 
to give, if anything, a larger value fordE/dT. That 
the reverse is the case would support the idea that 
the difference is experimental 
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To obtain double assurance on this point, however, 
a 4.98-mg radium preparation in a platinum-iridium 
needle 15 mm in length was measured both in the 
old [5] and in the new radiation balances. It was 
then encased in a glass tube of the same dimensions 
as the Hénigsechmid standards and remeasured in the 
new balance. In the first measurements the plati- 
num-iridium needle was entirely immersed in silicone 
oil in the cups of the old and the new balances. In 





the glass tube it had air around and above it, so that 
the temperature distribution should have approxi- 
mated that pertaining for the HOnigschmid standards 
The platinum-iridium needle used was radium source 
number 35917 referred to in Research Paper 2486 


; 

The value for the energy absorbed from radium 
source 35917 in the old radiation balance, corrected 
for the growth of polonium-210, was 747.5) pw, 
compared with 746.3 aw previously obtained [5] 
The values in the new balance, both unsheathed and 
sheathed in a 38-mm long glass tube of 3.1-mm 
internal diameter, and corrected for the growth of 
polonium-210 and to the equivalent in gamma-ray 
absorption of the old balance, were, respectively, 
749.8 uw and 749.1 ww The difference between 


these two values is insignificant 


5. Experimental Results 


During the first week of February 1954 the radia- 
tion balance was used to determine the rate of energ’ 
production of four national radium standards. The 
results of these measurements have already been 
reported in detail 3] In addition to these results 
it is, however, possible to derive a new value of the 
rate of energy production per milligram of radium 
element using the masses of the three standards that 
were determined by Hémgschmid on June 2, 1934 
The mass of radium element in the Canadian stand- 
ard is only a derived mass based on comparisons with 
the 1911 Paris and Vienna standards, and it is, more- 
over, doubtful if anv absorption correction was made 
in those comparisons. Only the masses determined 
by direct wewhing by Hoénigschmid have therefore 
been used in a redetermination of this constant for 
radium 

The masses of radium clement in the three primary 
standards as determined by Hoénigschmid on June 
2, 1934, and also the same masses corrected to 
February 2, 1954, using a half-life of 1,620 vears, 
are shown in table 3, together with the rate of energy 
absorption as determined by the radiation balance 
and the ratio of energy absorption in microwatts to 
mass of radium element in milligrams on February 
2, 1954 , 

The average value of the ratio of the rate of energy 
absorbed in the radiation balance to the mass of 
radium element, obtained by dividing the sum of all 
three rates of energy absorption by the sum of all 
three masses is 165.83 uwime of radium element, 
the standard deviation of the individual results being 
O18 percent 

Using the expression derived by Curie and 
Yovanoviteh [15], and later modified by Sanielevici 
[16] and Jordan [5], the growth of polonium-210 
in 19% vears, from June 1934 to February 1954, 
is calculated to give mse to a rate of energy produc- 
tron equal to 13.69 uww/mg of radium element The 
corresponding energy increment due to radium E 
is ‘0.87 uw/me of radium element. Hence the ratio 
of rate of energy production to mass becomes 151.36 
pw/me of radium element, when corrected "for the 





Tarte 3 Results for three national primary rad 
standard 


british U.8 U.s 
primary pramary prima 
radium radium radiu 
standard standard standa 

5432 5440, 43 

XV) XI\ 
Mass in milli 
grams, of radium 
element June 2 
1934 15. 60 20. 45 38. 2 


Mass, in milli 

grams, of radium 

element Febru 

ary 2, 1054 15S. 46, 20. 27, 37.4 
Rate of energy ab 

sorption in 

microwatts 2500.8 4300. 7 6203 
Microw atts per 

milligram of ra 


dium element 166. O1 165. 60 165. 80 


growth of polonium-210 and of radium E, or 130 
eal g-' hero! (all the polonium-210 alpha particles ; 
radium E beta particles being absorbed 

This result was, however, obtained with COL 
cups and glass tubes that together have a vamima-ray 
absorption equivalent to only O.184 em of lead 
Approximately 7 pereent of the energy dissipation of 
radium occurs in the form of gamma radiation. th 
absorption of which has been measured calor 
metrically by Zlotowski [17] up to an equivalent 
of 8 em of lead, at which thickness of absorber th: 
curve is essentially flat [5] 

Using Zlotowski's absorption eurve the value 
of 130.20 eal g°' hr-' for absorber equivalent to 
0.184 em of lead ean be corrected to “infinit 
thickness of absorber: i. € to an absorption quis 
alent to 8 em of lead A value of 138.6 eal eo! hy 
is then obtained 

The value of 130.20 eal go! hr 
have an accuracy of +0.5 percent 
deviation of the individual results is, however, as 
low as 0.13 pereent. Sources of radiation can 
therefore be compared with each other with a greater 
accuracy than their absolute energy emission can 
be measured, on account of the greater uncertainty 
in the determination of the effective value of dE/a7 


is estimated to 
The standard 


The old radiation balance [5] gave a value for this 
ratio of 128.9 cal ge 'hr-'! for 0.126 em of lead equir 
alent,’ which becomes 137.5 cal go! hero! when 
corrected to an absorption equivalent to 8 em of 
lead. This value was obtained, however, using 
needles whose radium content had been determined 
by gamma-ray comparison with the United States 
primary radium standard and not directly by 
weighing. The result obtained by Zlotowski [17 
was 139.6 cal g-' hr-' for an absorption equivalent 


to 8 em of lead. 





' The value quoted [5] was 128.9 calle~' br~' but this becomes 128.8 eal 5 
when the small radium E correction ts included 
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four national radium standards [3] 
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volt at 25° ¢ 
and 0.09 millivolt at 95° C 
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were 


Phe result 


trons of 


1. Introduction 


The silver silver chloride electrode Is employed 
extensively in the determination of ionization con 
tants and other thermodynamic data by the electro 
motive-force method [1] It is therefore important 
that the standard potential of this electrode be 
known as accurately as possible over a wide range 
of temperature 

Klectromotive-force measurements of cell A 





Pt: H, (2. 1 atm), HCL Om AvCl: Ag, (A 
at values of m suflicrently low to be useful in deter 
mining the standard potential by extrapolation to 
molality have made by a number 
investigators [2 to 16] The measurements of Gin 
telberg were made at 20° CC, and all of the other 
except that of Harned and Ehlers 


which covered the range 0° to 60° C, were confined 


zero been ol 


investigations, 


to 25° C. Re onthy, Harned and Paxton [17] have 
calculated the standard potential for the range 0 
to 50° © from the electromotive foree of cells of 


type A containing aqueous mixtures of hydrochloric 
acid and strontium chloride. In connection with 
the establishment of pH standards, the standard 
potential was needed in the range 60° to 95°C. In 
view of the extensive use of this electrode in electro 
chemical studies, it was deemed desirable to redeter 
mine the standard potential at lower temperatures 
as well a aad 

The measurements reported here were made at 17 
temperatures from 0° to 95° C and were limited to 
molalities between 0.001 and 0.12. The number of 





Figures in brackets indicate the literature references at the end of this paper 

The calculation of the standard potential from the data of Harned and Ehlers 
has been examined by Harned and Wright [10], Prentiss and Scatchard [11), 
Hamer, Burton, and Acree [12], Hills and Ives (14), and Swinehart [14] 
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from 
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Standard Potential of the Silver-Silver-Chloride Electrode 
from O° to 95° C and the Thermodynamic Properties 
of Dilute Hydrochloric Acid Solutions 


of the 


have 


and 
acid, and 


method of least 
were 
tandard potential wa 
0.02 millivolt 
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1) the 


Roger G. Bates and Vincent E. Bower 


eell without liquid junetion 
AeCl: Ag 
the 


been made of (1 tandard potential of 


2) the activity coeflicient of hydrochloric acid in aqueou 


from OF to 


relative 


oo ©, (3) the relative partial 


partial molal heat capacity of 
quares with the aid of punch-eard 
temperature, and SI cell 
found to be 0.22234 absolute 
at O° CC, O.OL millivolt at 25° ©, 
are compared with earlier determina 
derived from the 


analyzed at each 


to HOF ¢ 


electromotive force, 


cells studied ranged from 24 at 45° C and 55° C to 
SO at 60° Co and St at 25° CL The equations used 
for extrapolation were obtained by the method of 
least Puncheard techniques aided in the 


calculation 


Squares 


2. Experimental Procedures 


Hydrochloric acid of reagent vrade Wiis distilled 
in an all-glass still; the middle fraction (about two 
thirds) of the distillate was collected and redistilled 
The middle fraction of the distillate from the second 
distillation was diluted, as needed, with water to 
about O.l m and was standardized gravimetrically 
by weighing silver chloride, Test of the undiluted 
acid revealed no bromide [IS]. One of the three 
0.l-m stock solutions was standardized three times 
over a period of S months; the concentration ap 
peared to have changed only 0.02 percent in that 
time 

The cell solutions were prepared its needed by 
diluting portions of the stock solutions with water 


that had a conductivity of about 0.8% 107° ohm 
! 


1 


em at room temperature Dissolved air was re 
moved from most of the solutions by bubbling 
nitrogen; the rest of the solutions were saturated 


with hydrogen or boiled under vacuum. When the 
latter procedure was used, the weight of the solution 
was determined after boiling so that the final con 
centration could be calculated accurately. The 
electrolytic hydrogen, obtained in’ cylinders, was 
purified by passage over a platinum catalyst at 
room temperature and then over copper at 500° © 

Kach of the cells, described elsewhere [19], con 
tained two hydrogen electrodes and two. silver 
silver-chloride electrodes The latter were of the 
thermal-electrolytic type [2, 20]. The silver oxide 


from which they were prepared was washed 40 times 





with distilled water. The 1-\/ hydrochloric acid in 
which they were chloridized was a distilled sample 
free of bromide The electrodes were prepared at 
least 24 hours before use. For the high-temperature 
series (60° to 95° C the cells were provided with 
extra hvdrogen saturators consisting of three cham- 
bers, as deseribed by Bates and Pinching [21] 

Two calibrated potentiometers were used, The 
standards of electromotive force were a pair of satu- 
rated Weston cells maintained at a temperature 
near 36° C in a thermostated box of the type de- 
scribed by Mueller and Stimson [22 Three con- 
stant-temperature baths were emploved; water baths 
were used from 0° to 60° C and an oil bath from 60 
to 95° C. The temperature was regulated to the 
desired even temperature within the limits of +0.02 
deg C from 25° to 80° C and 0.03 deg C from 0 
to 20°C and above 80° C. Temperature measure- 
ments were made with a platinum resistance ther- 
mometel The difference of temperature between 
the oil bath and the solution in a cell immersed in 
the bath than 0.1 deg C at 
90° C 

The cells from which the data for the range 0° to 
60° C were obtained were measured initially at 25° C 
The constant-temperature water thermostat was 
lowered to near 0° C overnight, and the measure- 
ments from 0° to 30° C were made on the second 
lav, followed on the third day by the measure- 
ments from 30° to 60°C. A final check of 34 of the 
cells was made at 25 _&. The uverage difference 
between initial and final values was 0.18 mv. The 
final value was almost alwavs lower than the initial 
value, and there was indication that a con- 
siderable time was required for equilibrium to be 
established after the rapid drop from the higher 
Seven of the cells were measured 
only in the range 25° to 60° C. The data for the 
high range, 60° to 95° C, were obtained from a 
separate group of cells immersed in an oil bath. 
The initial measurements of these cells were made 
at 25° C or at 60° C, and the other temperatures 
were studied in ascending order. <A final check at 
60° C was sometimes but not alwavs made 

The electromotive-force values were corrected to a 
partial pressure of hydrogen of 1 atm. Inasmuch 
as the ionic strength did not exceed 0.113, the vapor 
pressure of each solution from 0° to 70° C was taken 
to be that of pure water [23]. The error introduced 
by this approximation appears to be less than 0.02 
mv at 70° © for the most concentrated solution 
studied. At 80°, 96°, and 95° C, the pressure cor- 
rection was made with sufficient accuracy by as- 
suming that the relative vapor-pressure lowering due 
to the presence of hydrochloric acid is the same as at 
25° C [24]. 

Hills and Ives [25] have identified an excess pres- 
sure effect due to the depth of the jet through which 
the hydrogen enters the solution. From their re- 
sults, it is evident that the effective partial pressure 
of hydrogen at an electrode located just below the 
surface is greater than that in the gas phase by (0.4 
A/13.6) mm, where A is the depth in millimeters of 
the hydrogen jet below the surface. In the cells 


was found to be less 


some 


temperature 





40 mm. 
o> Bed 


used in this work, h was about 


ah 
Che 


rection therefore amounts to 0.02 mv at 25° C. 9 


my at 6€° C, 0.08 mv at 90° C, and 0.16 my at 95 
Nevertheless, the corrections were not applied to 

electromotive-force data and standard potent 
reported here, in order that these results could 
used directly in other studies where the average 

depth is about the same (namely, 4 em) as in 

investigation. The thermodynamic constants 

hydrochloric acid solutions are unaffected, as tl 
depend upon the difference —F° and its cha 
with temperature. 


3. Standard Potential of the Cell 


From the equation for the electromotive forces 
of cell A one can write 


(log Mme log Y 


$H051SRT 
a F 


where /° is the standard potential of the cell, y~ is 
the stoichiometric mean ionic molal activity co 
efficient of hydrochloric acid, and the other symbols 
have their usual significance. Harned and Owen 
(1, chap. 11] have shown that experimental activity 
coefficients of uni-univalent strong electrolytes up 
to 1 m can be expressed with high accuracy by an 
equation of the form 


Aye 
log 7 \ -( ¢ > eat 


1+ Ba* ce 


log (1+-0.03604 m 


where ¢ is the molar concentration, A and B ar 
constants of the Debve-Hiickel theory, (’ is an ad- 
justable parameter, and a* is the ion-size parameter 
and (ert.) represents the total contribution of the 
extended terms in the Debye-Hiickel theory. 


When m does not exceed 0.1, ye differs from 
ymd°, where d° is the density of pure water, by 
less than 1 part in 1,000. Substitution of md 
for ¢ in eq (2) and combination with eq (1) gives 

A’ym 


E°"' = E°—gm Fant. | log m- 


t-(ext.)—log (1+-0.03604 m '} 


1+ B’a*yn 


where 8 is a constant for a particular temperature 
and value of a*. The values of A’ and B’ from 0 
to 100° C have been tabulated elsewhere [26) 
and (ert.) from 0° to 60° C for a*=4.3 is given by 
Harned and Ehlers [9]. The latter is only —0.00075 
at 0° and —0.00094 at 60° for the highest concen- 
tration studied in this investigation; hence, its 
value for 70°, 80°, 90°, and 95° C was obtained by 
linear extrapolation. These values of (ezt.) were 
used in the calculations at all the temperatures 
studied. The extended terms correction becomes () 
at m=0, but is a function of a*. The differences, 
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ert. (4.3 A ext. (6.0 A), at 60° C (where the best 
fit was obtained with a*=6.0) were not quite linear 
with m. Nevertheless, the mean departure from a 
straight line was less than 0.03 mv, or about 
one-third the probable error at this temperature 
The values of 2 30259RT I in absolute valts were 
computed from R=8.31439 j deg mole and 
F—96493.1 coulombs equivalent [27], and the 
absolute temperature, 7, was taken to be ¢ C- 
273.160 

The number of selutions studied was sufficiently 
large to justify the use of statistical procedures in 
analy zing the data With the proper choice of a”. 
a plot of £°’’, eq (3 hould be a straight line with 
intercept £° and slope — 8. The best value of a* 
is presumably the one that makes /°’’ most nearly 
a lmear function of » lo ascertain this best 
value, °°’ was calculated for three values of a* 
at O°, 25°, and 60° C and fitted to a linear equation 
by the method of least squares The standard devi- 
ation, ¢, of an experimental point from the least 
square line is plotted as a function of a* in figure | 
The curves are belheved to justify the selection 
of 4.3 A for a* at 0° and 25° C and 6.0 A at 60° C 
The values of a* tor temperatures between 25 
and 60° and from 70° to 95° C were determined by 
inspection of the plots of /°’ as a function of m 
for two or more values of a* 

If an incorrect value of the ion-size parameter is 
ised, the plots of i = with respect to m become 
curved, and the intercept of the straight line estab- 
lished by least squares is no longer the true value 
of bk The influence of a change in a* is demon- 
strated by the data for 25° C 


1 E 7 
\ mit 
2.0 0. 22222 0.19 
1.3 22234 07 
6.0 22246 13 


Table 1 contams a summary of the least-square 
calculations at the 17 temperatures. The standard 
potential of cell A is given in the fifth column. 
The standard deviation, ¢,;, in millivolts, of the 
intercept is given in the sixth column. The value 
of £° from 0° to 90° C is given by the equation 


ip 0.23659— (4.8564 * 107") t—(3.4205 * 10~°) t? 


+-(5.869 > 107%) t, (4) 


where fis in degrees Celsius. The standard potential 
of the silver-silver-chloride electrode is either equal 
to £° (cell A) or —E°, depending on which of the 
two common conventions for single electrode po- 
tentials is adopted. 

The “observed” values of E° are compared in 
table 1 with those calculated by eq (4). The lasi 
column gives A, the difference in millivolts, between 
the calculated and observed value at each tempera- 
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ture. The average value of A at the 16 tempera 
tures is 0.04 my. 

Figure 2 is a plot of £°’’ at 0°, 25°, 60°, and 90° C 
(open circles) as a function of molality. The closed 
circles were computed from the data of Harned and 
Ehlers [9] by the method described above. They 
lead to values of 0.23660 abs v for E° at 0°, 0.22252 
v at 25°, and 0.19650 v at 60°. 
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4. Activity Coefficient of Hydrochloric Acid 





Che electromotive forces given in table 2 were com- 
puted from the smoothed values of /°”’ at round 
values of the molality This # was computed, 
in turn, from the intercepts and slopes of the least 
square lines listed in table 1 The mean activity 
coeflicients calculated by eq (1) from these smoothed 
values of / and the values of /° given in table 1 are 
summarized in table 3 

Neither the electromotive force nor the activity 
coeflicient was smoothed with respect to temperature 
Henee, for a calculation of the thermodynamic quan 
tities derived from = the temperature coellicients of 
electromotive force, the values of 





log y, at 25°C 





and at intervals of 10 deg (* from 0° to 90° were 
fitted by the method of least squares to a power 
series in f, the temperature on the Celsius seale 

log y, = A+ Bt+Ct (5 


The values of log y, were given equal weight at each 


temperature. The constants of this equation for 
eight values of the molalitv are listed in table 4 
lanie 4 Constants of the equation log v4 1 it Ct 
for the te mperature range t i fot ww) ¢ 
! betw calculated and observed value 
n percent { lon ‘ wy { 
1 / f s 
Percent 
, aor 0 O1470 0. 273K" 0. 27™« Ww O30 
on C2051 ISA 1.30 24 
on O10 445 10 “4 
! O40 10 13 4 
" 4 ‘ san wv 
wee 1. O08 ‘ ‘ Ww 
Bl 1a ~ is 
re ] 11 “4 


The last column vives the mean difference between 
at the Ll tempera 

tures, expressed as percentage of log ¥, at 25° C 
' When the values of log YY, were weighted according 
to the reciprocal of the probable error of [2°"" at 
the appropriate temperature, the fit to eq (5) was 
not as complete as when equal weight was given to 
each value. The relative partial molal heat content 
computed from the two sets of constants differed on 
the average by 15 j mole! at 0° C, 7 j mole”! at 
25°C, and 42j mole"! at 90°C. The relative partial 
molal heat capacity was changed about 0.4 j) deg”! 
mole! at 0° C. 0.5 } deg™' mole™' at 25° C, and 1 2 j 
deg”! mole~! at 90° C 


the calculated and observed log v4 


, 


5. Relative Partial Molal Heat Content and 
Heat Capacity 


The temperature variation of log Y, can be used 


to caleulate the partial molal heat content, L,, and 


partial molal heat capacity, J,, of hydrochloric acid 
relative to its value in the iniinitely dilute solution. 








Closed cirel 
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The former is given by 
2) low ¥, Ls & 
“ —~, (ty) 
o7 1.6052R7 
where 7 is the temperature on the Kelvin scale 


Inasmuch as O7' Of, we obtain, by combination of 


eq (5 and (6 
L.A H052RT(B 4 204 (7) 
and 
J oes O.2104ART°C4 9 210ART( B+ 201 (S) 
o7 


The values of L, and .J,, 


from these two equations are listed in tables 5 and 6 


in absolute joules, calculated 


The relative partial molal heat content at 0 
25°, 60°, and 90° C is plotted as a function of m' 
in figure 3. The dots represent the results obtained 
by Harned and Ehlers |9, 1] at 0°, 25°, and 60° C 


The dashed line locates Sturtevant’s calorimetric 
values at 25° C [28]. The agreement with the 
earlier determinations can be regarded aus very 


satisfactory at O° and 25° C and acceptable at 60° C 


The relative partial molal heat capacity, J,, at 25° C 
is plotted in figure 4. The dots again indicate the 











values obtained from the measurement of Harned 
and Ehlers The dashed line is an extension to 
These are the means of the two sets of values given by Harned and Owen 
computed from the experimental data in two different w 4) 
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Fiaure 3 Relative partial molal heat content. L at hydro 
chloric acid at 0 Zi) 80°” 
aquare root of the molality 


and 90° C as a function of the 


of Harned and Eble Dashed line repre 
calorimetric results at 2 ‘ 


indicate the reeult 
ents Sturtevant 
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Discussion 


The values of /° for the temperature range 0 
60° C are compared in table 7 with those obtain 
from the measurements of Harned and coworke: 
9, 17). The standard potentials of Harned at 

. Paxton, given in the fifth colemn, are in’ bett 
eo T ; agreement with the present work than are those 
Harned and Owen (second column Although thi 
values are based on only six points below an jon 
strength of O.1, Harned and Paxton point out th: 
a straivht line could be drawn to within 0.08 my o 
these six points at nearly every temperature 
Harned and Wright's recalculation [LO] of Harn 
and Ehlers’ data, based on improved values of thy 
natural constants, lowered the figures in the second 
table 7 by an average of about 0.14 mm) 
whereas Swinehart’s recent reealeula 





column of 
5 ' O° to 40° C 








, 
: te 
tion [14] with the aid of newer values of 2, 7, and / 
raised them by 0.09 mv on the average The 
extrapolation method of Harned and Ehlers, which : 
“ expressed the activity coeflierent in eq (1) by the 
Debve-Hieckel limiting law, was used to obtain all of ) 
; the potentials except those given in the last column | 
In the present investigation it was found that such | 
o , an extrapolation procedure, applied to data at 
molalities up to O.1, vields a curved line, concav: 
a , upward, and of appreciable slope at low concentra 
° tions I 
" 
Tanue 7 Standard potentia / of ce | am O° ta? ( ( 
¥ n absolute volt I 
a 
EEE ——EE ——————ee 
7 5 o-4 | 
mire Electromotiveforce data of Harned at li 
Ehlers [v 
hn t ne ! Relative } un | ul heat pa j ofl ire ‘ — r? ( 
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Chwen [1 recalculated culated 
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‘ at ft il net lata obta ! oke ‘ hheminke n 
sbowe } ‘ 
0 oO. [a2 Oo. 285 0. 204 Oo. °ONN? 0. 2 nt 
Moo wwe PAD PAT LN) Mis 
i) Nis wi 24145 Jaa) 142 (; 
1S , TIM wth Js 22857 } 
2» YMA 2a v2 PIKS7 i 
nl 
- onnen o77e0 one eons 99034 
m' 0 of the straight line representing the values * oa oaea oe yee ase rn 
. , hee Saas 5 “ 207 Hise "216 21207 21208 
) rb . ) . ) , 22 2 
of J, obtained calorimetrically by Gu ker and . oneal en ono od — a 
Schminke [29] * at molalities from 0.1 to 2.25 
“ yea S2 JAM ay 0 
si AS 20042 200% 20045 it 
‘The data of Gueker and Schminke deviate sharply from the straight line wo 1820 10645 19640 ! 
helow O11 @ his anomalous behavior has net been exphined \\ 


















Evidently the consistency of the different sets of 
data ean only be judged if both sets are treated in 
the same mannet As may be seen in figure 2, the 
electromotive-foree data and standard potentials re 
ported here are im acceptable ayvreement with those 
of Harned and Ehlers at 0° and 60° C, but appear 
to be about O.1IS mv lower at 25° C A difference 
of this magnitude at 25° C, where the results are 
siatistically the most precise, is difficult’ to explain, 
particularly because the silver silver-chloride elec 

trodes and the hydrochloric acul were prepared by 
smnilar two investigations \ 
critical examination of the eleetromotive-force data 
obtained by other workers is therefore of particular 


procedures in the 


interest 

This comparison was made first at low concentra 
tions, Where the mode of extrapolation has the small 
est influence on the result All of the available emf 
data were accordingly converted to absolute volts 
by multiplying by 1.00033 [S80]. Values of °°" were 
with a@*--4.3. The results 
below OOO} ar 
The open circles are the data of 


then computed by eq (3 
of this reealeulation at > molalities 
shown uw figure 5 
this investigation, and the least square line is showy 
The dashed line is the extension of the straight line 
through the points of Harned and Ehlers, all of which 
were at molalities above 0.003 

The most numerous data im this region of low con 
Anderson and Young [15] 
The value 
measurements appears to 


centrations are those of 
madieated by closed circles in the figure 
of / 
be about 0.22242 abs v The crosses were calculated 
Carmody |S], and the 
half-shaded circles mark the lowest points of Roberts 
i The other data for cell A in’ this low-range 
display larger deviations and are not plotted, The 
four measurements of Linhart [8] m the range of the 
figure, all below 0.001 m, vary from 0.2225 to 0.2228 
The five pomts of Maronny and Valensi [16] below 
0.0025 m le O01 to 0.4 mv below the solid line The 
value of £° computed from Nonhebel’s six 
measurements [5] between m—O.0008 and m—0.008 
is 0.22943 4. 0.00005 abs v Below m— 0.0008. how 

ever, rises rapidly, exceeding 0.225 v at the 
lowest molalities studied 


a 
obtained from. these 


from the measurements of 


average 


\ comparison limited to low concentrations suffers 
from the fact that the experimental data are usually 
less accurate below 0.01 m than above Henee, the 
eleetromotive-force data of Giintelberg [6] at 20 
and of Roberts, Carmody, Harned and Ehlers, and 
\nderson and Young at 25° C for molalities up to 
Ol m were smoothed to round molalities, where 
necessary, on a plot of #° as a funetion of m and 
are compared in table 8 It is seen that the values of 
CGintelberg agree reasonably well with those reported 
here and are somewhat lower than those of Harned 
aml Ehlers. The latter are also higher than the 
others at 25° C, whereas those of Carmody and of 
\nderson and Young agree well with the present 
work The emf data obtained by Roberts between 
0.01 m and 0.1 m appear to fall between those of this 
mvestigation and the data of Harned and Ehlers 
With the exception of one low value, obviously 
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a, 


Fiaure 5 


erroneous, the live pots of Noves and Ellis || below 
m O10 agree with the results reported here, as do 
the four of Seatehard [4] between 0.0104 and O01 m 
and the two of Linhart [3] above 0.01 mm. Seatehard’s 
three points near 0.01 m, however, lie nearly 0.2 my 
below the line through his other points. It may be 
concluded that the work of Giintelberg, Carmody 
and of Anderson and Young is consistent with the 
present study, whereas the measurements of Non 
hebel and Roberts, and those of Linhart below 0.01 
m, tend to support the higher value of Tlarned and 
Ehlers at 25° C 
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\ rather large upward trend in & at the lowest 
concentrations is observed in the data of Linhart as 
well as of Nonhebel \ departure from the theoreti 
cal slope is not to be expected in this region and was 
not found by Carmody or by Anderson and Young 
It is possible that traces of oxygen, known to shift 
the potential of the silver-silver-chloride electrode 
toward more positive values in acid solutions, may 
explain this elevation of electromotive force at low 
molalities. The chloride-ion concentration in’ the 
vicinity of the silver-silver-chloride 
lowered by the following reaction [6] 


electrode is 


2 HCO 


2 AvCl+ LO (DY) 





The resulting change of emf may be appreciable in 
dilute solutions, for dE /dn.,, where ney is a number 
of equivalents of chloride ion, is much larger than in 
solutions of moderate or high concentration 

Nevertheless, dissolved air cannot explain the dif 
ference between the results of the present investiga- 
tion and those of Harned and Ehlers because an 
air-free technique was used in both investigations 
The potentials of silver silver-chloride electrodes are 
known to be altered likewise by traces of bromide 
118, St} and by aging during the first 30 hours after 
preparation [$2]. A lowering of the electromotive 
force of the cell by 0.18 mv would require about 
0.02 mole percent of bromide impurity in the hydro 
chloric acid used in this study; this quantity could 
hardly have gone undetected in the test that was 
performed. The effect due to aging causes the emf 
of the cell containing a freshly prepared. silver 
silver-chloride electrode to be too high. The agree- 
ment among measurements at 25° C made at dif- 
ferent points in the temperature series would seem 
to rule out a pronounced effect due to aging. No 
simple reasonable explanation for the differences 
between emf values at 25° C reported here and 
those of Harned and Ehlers has been found 

The activity coefficients and other thermodynamic 
properties of hydrochloric acid are dependent not 
upon the value of /° but on the difference /—F 
Inasmuch as the extrapolation lines are nearly par- 
allel (see fig. 2),.° the activity coefficients at 25° C 
reported here agree very well with those computed 
from Harned and Ehlers’s measurements with a 
standard potential of 0.22252 abs v (the value ob- 
tained from the emf data of Harned and Ehlers by 
the extrapolation procedure used in the present 
work The activity coefficients from these two 
sources are compared in table 9 with those obtained 
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by Hills and Ives [13] in a careful study of the 
hvdrogen-calomel cell without liquid junetion and 
with those computed by Shedlovsky [33] ° from trans- 
ference numbers and the electromotive force of cells 
with transference. The agreement with the deter- 
mination of Hills and Ives is very satisfactory, and 
the only notable difference from the values of Shed- 

’ Prentiss and Seatchard [11] have noted that the slopes of the lines plotted 


from the data of Carmody, Roberts, and Harned and Ehlers are nearly the same 
* A similar comparison has been made by King [4 


lovsky appears to be at m 
corresponds to 0.19 mv in the electromotive fo 


0.1, where the depart 
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for advice and assistance in the statistical treatm 
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are also indebted to T. F. Young for helpful dis: 
sions during the course of the work 
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The phase diagram of the gold-uranium 
thermal analysis 


The Gold-Uranium System ' 
R. W. Buzzard and J. J. Park 


metallographic examination, and X-ray 


svstem Was constructed from data obtained by 
diffraction The system i 


characterized by two intermetallic compounds one forming peritectically at 1.216° C and 


having an apparent composition of U,Au 
mately 1,450° C 
one between U Au and gold at 855 C and &7 


and having an apparent composition of UAu 


and the other melting congruently at approxi 
There are two eutecties 
S atomic percent of gold, and the other between 


uranium and U,Auy at 1,105° C and 10.5 atomie percent of gold The solubility of uranium 


in gold appeared to be about 0.6 atomic percent at 855° C, and the 


uranium Was approximately 3.2 atomic peres 


transformation temperature of uranium from 762° to 738° C 


tion from 653° to 647° ¢ 


1. Introduction 


lhe purpose of this investigation was to make a 
survey of the gold-uranium system and to develop 
the constitution diagram by correlation of thermal, 
X-ray, and microscopic data 


2. Previous Work 


\ survey of the literature ? showed that’ prior 
information on the constitution of the gold-uranium 
svstem meagel Metallographic studies of 
alloys having compositions of 5 to 12 atomic percent 
of gold have been reported.4* The alloy contain 
ing 5 atomic percent of gold, in the “as-cast’’ condi 
tion, showed a definite two-phase structure, and the 
alloy having the composition of 12 atomic percent of 
vold, in the as-cast condition, showed particles of a 
gold-rich phase in a ettectic matrix, [t was reported 
that some solubility of gold in alpha uranium was 
indicated 


was 


3. Preparation of Alloys 


The component metals consisted of uranium of 
about 99.9-pereent purity (Mallinckrodt Biscuit) and 
W0-fine gold.° Most of the alloys were prepared in a 
high-frequency induction furnace under vacuum, 
using either beryllia or thoria crucibles.’ The cold 
furnace was pumped to about 5 u, and on melting, 
the pressure attained was about 200 u. The alloys 
in the composition range 30 to 99 atomic percent of 
gold reacted with the beryllia crucibles to the extent 
that the beryllia adhered to the ingot, whereas the 
alloys in the range 0 to 30 atomic percent of gold slid 
free of the beryllia crucibles, giving no visual indica- 
tion of a reaction of the alloy with the refractory 
Subsequent metallographic examination (fig. 1, A) 
revealed that the alloys in the range 50 to 70 atomic 
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solubility of gold in 
Gold lowered the gamma-beta 
and the beta-alpha transforma 


nt at 1.105° C 


old were partially oxidized, apparently 
from a reaction with the beryllia crucibles. This was 
later confirmed by observations during thermal 
analysis of the alloys in this composition range 
These alloys were completely oxidized when bervilia 
crucibles were used, the resultant oxide mass con 
taining pellets of practically pure gold. The use of 
thoria crucibles for the melting and thermal analysis 
of the alloys in this range of composition apparently 
avoided the oxidation of the alloys, although it was 
evident, on visual examination of the alloys, that a 
reaction between the alloy and the refractory had 
occurred, Metallographic examination of the thoria 
melted alloys showed no evidence of the oxidation 
that was observed in the beryllia-melted alloys, but 
chemical analysis revealed that the thoria-melted 
alloys picked up as much as 1.6 weight percent of 
thorium. Subsequent melts of alloys in this com 
position range were made by are-melting 

The exploratory series of alloys were made from a 
50 weight-percent master alloy that, on microscopic 
examination, appeared free of oxide. Subsequent 
alloys were made by direct alloying. Chemical 
analyses were made for uranium on the as-cast alloys 
and the gold determined by difference. When visual 
evidence of attack on the crucible was observed or 
microexamination gave evidence of oxidation, analy- 
ses were made for both component metals. The 
compositions of the alloys used are given in table 1. 
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Fieurs 1. Gold-uranium alloys 


x 1 


A. (5 atomic percent of gold alloy, inchoate oxidation (grey), UsAus (light), unetched 
B. 89.8 atomic percent of gold alloy; eutectic-ty pe structure of gold phase (light) and epsilon (dark 
eutectic structure, containing excess delta (light) outlining grains, unetched “ 

swcetic acid electrolytic etch 


5 percent of KC N electrolytic eteh 
C. 14.5 atomic percent of gold alloy x 500 
LD). 65 atomic percent of gold, epsilon phase (dark) outlining grains of delta, orthophosphoric acid plus 
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4. Procedures 


4.1. Thermal Analysis 


Thermal arrests were determined by time-temper- 
ature analysis, ina molybdenum-wound furnace. In 
eases where an uncertainty of a thermal arrest 
existed, it was checked by differential analysis 
The accessory apparatus included an air-operated 
controller to maintain uniform rates of heating and 
cooling and an electronic recorder to reproduce auto- 
matically the heating and cooling curves The 
thermal curves were obtained from = specimens of 
approximately 50 g, which were heated and cooled 
at controlled rates of 2 deg C/min, with the furnace 
operating under a vacuum of approximately LO 4 
The thermal arrests were derived from the cooling 
eurves 


4.2. Microscopic Analysis 


The specimens for microscopic analysis were 
mounted in Bakelite, ground on a series of silicon 
carbide papers, with the finishing paper of 000 grit, 
and finished with levigated alumina either on a wet 
broadcloth lap or a dry paraffin lap. Excessive 
oxidation was encountered during the hot pressing 
of the Bakelite and on subsequent wet polishing of 
the allovs in the range 60 to 75 atomic percent of 
gokl. These alloys were mounted in dental acrylic 
casting resin and dry polished 

In general, the structures of the uranium-rich 
allovs were developed by electrolytic etching in a 
chromic-acetic acid solution (100 ml of glacial acetic 
acid, 20 g of chromic acid, and 30 ml of water) 
by 2 to 4 immersions of S-sec duration at 70 v and 10 
amp/em’, The structures of the gold-rich alloys 
were developed electrolytically in a 5-perecent 
solution of potassium cyanide by 2 to 4 immersions 
of 30 sec each at 4.5 v and 1 amp/em*®. The strue- 
tures of the alloys in the range of 45 to 75 atomic 
percent of gold were developed electrolytically in a 
solution containing 50 ml of orthophosphoric acid, 
100 ml of concentrated sulfuric acid, and 100 ml of 
water by 2 to 4 immersions of 15-sec duration at 
1Ovand 10 amp/em*. This latter etching technique 
subsequently produced equally satisfactory results 
on alloys of all compositions and was used exclusively 
in the later stages of the investigation. 

Specimens to be quenched were sealed under 
vacuum in 96-percent  silica-glass tubes of 7-mm 
hore, homogenized for S days, and furnace-cooled 
\ homogenizing temperature of 825° C was used for 
vold-rich alloys, and 1,050° C for the uranium-rich 
allovs. The homogenized alloys were subsequently 
reheated and quenched in ice water. The nickel- 
block technique was used in the manner previously 
deseribed ’ 


4.3. X-Ray Analysis 


The specimens used for the microscopic studies 
were subsequently placed in an X-ray spectrometer, 
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and the charts were obtained of the (Cu-Ka, or 
Co-Ka) X-ray diffraction maxima at room tempera- 
ture. Powder patterns were used to supplement the 
spectrometer data. By this method it was possible 
to identify the phases present in the alloys and to 
approximate the alloy compositions at which the 
compounds appeared. ‘These results were correlated 
with the thermal and microscopic data 


5. Results 
5.1. Thermal Data 


The results of the thermal analyses are given in 
table 2 and are also incorporated m the equilibrium 
chagram of figures 2 and 8. Additions of gold to 
uranium lower the melting point of uranium from 
1.183° to 1,105° C at 10.5 atomic percent of gold, 
from which the liquidus temperature rises to a 
maximum in excess of 1,450° C at 75 atomie percent 
of gold; from the maximum the liquidus temperatures 
fall sharply to 855° C at 87.5 atomic pereent of gold, 
and then rise sharply to the melting point of gold at 
1,060° C 

The results of the liquidus determinations between 
30 to SO atomie percent of gold were subject to an 
error due to a reaction between the alloy and refrae- 
torv. In this range of composition the allovs were 
very reactive in the liquid state, and during thermal 
analysis a reaction was noted between the alloys and 
the bervilia crucibles. Sometimes a reaction was 
observed when using thoria crucibles. With bervllia 
crucibles, the allovs showed a decided tendency to 
oxidize completely during the thermal-analysis cvele 


Prarie 2 Summary of thermal data obtained for gold-uranium 
alloys 
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the oxidized mass apparently consisted of uranium 
oxide that was interspersed with pellets of practically 
pure gold. The thermal reported were 
obtained from alloys tested im thoria crucibles and 
represent the arrests material that 
appeared, on subsequent metallographic examina- 
tion, to be of homogeneous oxide-free metal. The 
alloy having a composition of 76.1 atomic percent of 
gold apparently did not melt at 1,450° C; this was 
the limiting maximum temperature of the thermal- 


arrests 


observed on 


analysis furnace 


Liquidus arrests reported aby 
1300° C are approximations, as the characteristi: 
the furnace did not permit accurate determinati: 
in excess of 1.300° (C 

Thermal analysis of the alloys in the composit 
range 20 to 60 atomic percent of gold showed 
in the viemity of 1216° CC, wh 
indicated a temperature horizontal in this rang: 
composition, A second temperature horizontal \ 
observed in the range of 3.7 to 60 atomic percent 
gold in the viemity of 1,105° C 

The gamma-beta transformation temperature 
uranium was lowered by gold from 762° to 738° ( 
This value was constant for alloys containing 1 to | 
atemic percent of gold. The beta-alpha transform 
tion temperature was similarly lowered from 653° | 
to the constant value of 647° C' for the same allo 

In the gold-rich alloys (75 to 99 atomic percent 
wold) a 855° C Mh 
together with the lowermg of the melting pom 
gold to a minimum at 855° C peres 
of gold by additions of uranium, was indicative of 


consistent arrest 


consistent arrest occurred at 


at S7.5 atom 


eutectic transformation occurring in this cOompositio 
range 

No solid-phase reactions were noted below 1 216° ( 
in the composition range 60 to 75 atomic percent of 
gold, but it should be noted that only one allov wa 
examined in this composition range No indieation 
of the euteetic arrest was the allo 
contamimng 1.0, 2.5, and 99.7 atomic percent of gold 


obser \ ed inh 


5.2. Microstructures of the Alloys 


In the study of the microstructures of the gold 
uranium alloys a eutectic structure was observed in 
the range 75 to 99 atomic percent of gold (fig. 1, B 
In the range of 75 to 87.5 atomic percent of gold a 
second phase appeared in increasing proportions a 
successive increments of uranium were added to thy 
allovs. This second phase was optically 
under polarized light, and therefore it was concluded 
that this phase was not uranium. The alloy becany 
increasingly brittle as the uranium content 
raised an indication of an intermetallic compound 
Similarly, in the uranium-rich alloys in the range 2.5 
to 60 atomic percent of gold two phases were ob 
served, One phase consisted of a typical eutecty 
structure (fig. 1, C The second phase was optical) 
enisotropic and therefore differed from the second 
phase observed in the gold-rich alloys. This was 
indicative that at least two intermetallic compounds 
existed in the system 


iso tropr 


Was 


A two-phase structure (fig. 1, D) was observed in 
the range 60 to 75 atomic percent of gold. One 
phase was identified as that common to the gold 
rich alloys, and the second as that common to the 
uranium-rich alloys. From this evidence tt was 
concluded that apparently only two intermetally 
compounds exist in the gold-uranium system 

The identifiable phases observed in the micro 
structures of quenched alloys (table 3) were used in 
the determination of the uranium and gold solvus 
The maximum solubility of uranium in gold was found 
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Panne 3. Microscopie analysis of gold-uranium alloy the diffraction peaks of uranium were determined by 
lid solution; 2, two-phase structures; F, fusion comparison with known alpha-uranium peaks.” ‘The 
remaining unassigned peaks were attributed to a 
ear ee eae ae eee second phase (delta); the intensities and number of 
—_ { gold the delta peaks increased as the percentage of gold 
pera increased, whereas the intensities and number of the 
se ' - _ alpha-uranium peaks gradually decreased. How 
ever, the alloy containing 65.5 atomic percent of 
scandens gold gave a diffraction pattern consisting of the peaks 
; attributed to the delta phase plus weaker peaks later 
1, 0 ' considered to be from the epsilon phase. Similarly, 
a aoe , , in the region 75 to 99 atomic pereent of gold, the 
( nee ; ; ; strength of the known gold peaks decreased with 
increasing amounts of uranium, though no peaks 
br attributable to uranium itself were detected; the un 
( = known peaks of this series of alloys were considered 
to be from a second phase (epsilon) The pattern 
RIE ee Se et eee from the alloy contamig 76.1 atomic pereent of vold 
Quenel ; { gold consisted almost entirely of epsilon peaks. By such 
pera a comparison of intensities of known and unknown 
, — m6 | m2) om peaks, the composition of the delta phase was esti 
yf ; mated to be 60) atomic pereent of vold with the 
rT a formula U,Auy, and the pure epsilon phase was esti 
mated to consist of 75 atomic percent of vold with 
( 1, an ' the formula UAu 
tof —s The known d-spacings for gold at room tempera 
wit ; wa , ' ture were not changed appreciably by additions of 
Pion } 
lo : Ni Tanue 4 /nterplanar spacings of compounds * 
old : ROD 
"0 
6on Uranium Dette cpalien ch ae Epstion | Epsiion at 76 — 
plus delta spacing itomic percent pacing atomic percwnt | id at 
At Me atom most of gold most of gold nh 
old Acearding to thermal data these structures should be two phase; undoubted! percent of | probatle probable atone le 
doin econd phase was not developed by the metallographic technique used one wane Powder olid - Powder old no Sold. 
Ib 
ld it ‘ - , i i 1 1 ! 1 
to be in the range 99.2 to 99.7 atomic percent of gold +: a 
a The solubility of gold in alpha and beta uranium 7 yar eee Lo | 2 | 2M a 
was found to lie between 1.00 to 2.5 atomic percent 4 » ant “ ‘nt 
<8 of gold. The microstructures indicated that) the a _— we | ae 23 oe oy | 3 
ule solubility of gold in gamma uranium was in the ' 2 whe 
—_ vicinity of 3.7 atomic percent of gold. The plot of hom oe fee ee cee 
= &F the thermal-analysis data, however, located this at shen Isheedll indeed bt 
und. ' +2 atomic percent of gold | 
e 2.0 | t » wit 
rob i = ¢ 74 ; 7 712 ; 716 17 17 17 
wn 5.3. X-ray Data ihm | GAM | DAM ast | aan | 
5 al 
cally The X-ray diffraction patterns were obtained from nae 
cond the same surfaces that were used for the metallo ee 1a 
wes graphic examinations, and the data obtained were = ava | 147 
uns correlated with the microstructures. The alloys +. 08 1. 382 
quenched from the gamma field were identified by 14 1308 | 8001 OM ee ee 
1 usa 1 JAA 1 wa 1 2a 
ed in the alpha-uranium lines (transformed gamma) All \ asi 
One the alloys in the range 1.0 to 59.9 atomic percent of “me tome | cm lem ic Pe 
gold vold gave diffraction patterns characteristic of the 
» the alpha uranium and a second phase (delta); in the va |e an pleats 
wis range 65.0 to 74.1 atomic percent of gold a typical a. O58 oe | oe | oe 1‘ 
tall two-phase pattern was observed characteristic of OMS OO | 00M 0 fs 
delta and a second phase (epsilon); and the alloys oss jae ons 
cro containing in excess of 76 atomic percent of gold » wis 
ed in gave diffraction patterns typical of epsilon and gold 
yvus In the region 1 to 60 atomic percent of gold, the * Caloulated from the 29 value Uranium Cold 
found X-ray ditisactian chart were endipandpnd comoased: +C. W, Jaco and B, E, Warren, J. Am, Chem, Soe, 88, 2586 (1997 
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uranium, nor were the d-spacings of alpha uranium 
changed by additions of gold No appreciable solid 
solubility of either compound was detected The 
interplanar spacings of the compounds (table 4) 
were determined, but the crystal structures of the 
compounds were not ascertained 

Powder methods of X-ray analysis were not com- 
pletely successful owing to the high affinity of the 
compounds for oxygen. The microstructure of the 
alloy having a composition of 76.1 atomic percent of 
gold consisted almost entirely of the epsilon phase; 
the sample was so brittle that minute pieces (approx- 
imating a powder sample) were scraped off the ingot, 
sieved, and used as a powder sample. The many 
weak peaks obtained indicated a compound of low 
symmetry 


6. Proposed Diagram 


The thermal arrests, derived from the cooling 
curves, were obtained from as-cast alloys that were 
cooled rapidly from the liquid condition. Chemical 
analysis showed these ingots to be homogeneous 
from top to bottom. ‘Thermal arrests were observed 
at about 1.216°, 1,105°, 738°, and 647° C in most 
uranium-rich alloys (0 to 60 atomic pereent of gold) 
and at about 855° C in most gold-rich alloys (75 to 
100 atomic pereent of gold); in the range of 60 to 75 
atomic percent of wold, no thermal arrests were noted 
below 1,216° CC. Metallographic studies revealed a 
eutectic-like structure in both the uranium-rich and 
gold-rich alloys and a two-phase structure in the 
range 60 to 75 atomic percent of gold Both X ray 
and metallographic analysis indicated the existence 
of two compounds having apparent compositions of 
U,Au, (delta) and UAuy (epsilon The reaction 
horizontal observed by thermal analysis at 1,216° C 
is peritectic in nature, denoting the formation of the 
delta phase by this reaction 

In the gold-rich alloys, correlation of the metal 
thermal analysis data indicated a 
eutectic at 855° C and approximately 87.5 atomic 
percent of gold. Similarly, in the range 4 to 60 
atomic pereent of gold, correlation of metallographic 


lographie and 


and ‘thermal-analysis data indicated a eutecti: 


1,105° © and approximately 10.5 atomic perce: 
of gold. The solubility of gold in uranium yw 
established by correlation of metallographic a) 
thermal data at approximately 3.2 atomic perce: 
at 1,105° C, and the solubility of uranium in go 
was established at approximately 0.6 atomic perce 
at 855° C (figs. 2 and 3). Two compounds, havi: 
the apparent compositions of 60 and 75 atomic px 
cent of gold, were identified in the system. ‘T! 
compound at 60 atomic percent of gold was form, 
peritectically at 1,216° C, and the compound at 7 
atomic percent of gold) was congruent meltin 
with a melting point in excess of 1,450° C 


7. Summary 


hig. 3 


The gold-uranium system was constructed 
from data obtained by thermal, microscopic, and 
X-ray analyses 

The system is characterized by (a) two eutecti 
one occurring at 1,105° C and a composition of 10 
atomic percent of gold, and the second at 855° (© 
and 87.5 atomic percent of gold, and (b) two inte 
metallic compounds, one formed peritectically at 
1,.216° C having an apparent composition of 60 
atomic percent of gold, and the second melting con 
gruently in excess of 1,450° C, having an apparent 
composition of 75 atomic gold. The 
solubility of uranium in gold appeared to be between 
0.3 to 0.6 atomic percent at 855° C, and the solu 
bility of gold in uranium was approximately 
atomic percent at 1,105° CC. The gamma-beta trons 
formation temperature of uranium was lowered from 
762° to 738° C and the beta-alpha transformation 
temperature was lowered from 653° to 647° C by gok 


percent of 


The authors express their appreciation to Martha 
S. Richmond and J. R. Baldwin for the chemical 
analyses 
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A Study of Absolute Standards of Mutual Induct- 
ance and in Particular the Three-Section 
National Bureau of Standards Type 


Frederick W. Grover ' 


rhe results of a study of the number and location of the cireles of zero field surrounding a multiseetion 


coil are presented The configuration of the field surrounding the equatorial region of several three-seetion 
coils has been partially mapped out One arrangement vields a design in which the mutualinductance con 
tribution of anv secondary turn differs but little from that of the median turn By locating the secondary 


vmmetrically about the cireles of zero field, the correction for winding distribution ean be made small, and the 
effect of uncertainties in location of individual secondary turns minimized 


1. Introduction 


Certain absolute measurements, as for example, the absolute measurement of resistance 
depend upon the use of a mutual inductance whose value can be calculated from its geometrical 
dimensions and the number of turns on primary and secondary 

In order to minimize the effect of dimensional uncertainties the standard should be of 
fair size with relatively large separation between primary and secondary, Further, to obtain 
sufficient sensitiveness in the measurements in question, the value of the mutual induetance 
should be as great as of the order 10 millihenrys. These conditions demand a value of the 
product N.N, of primary and secondary turns of the order of 100,000 

The dimensions of a single lave coil, wound in a serew thread on a suitable evlindrn al 
form, may be measured with great precision. With such a large requisite value of NVA 
however, it is impracticable that both the primary and secondary should be single-laver coils 
Of necessity, one winding should be of the form of a multilaver coil, wound in a channel of 
square or rectangular cross section 

In a form of mutual-inductance standard designed by Albert Campbell in) 1907 [1]? and 
used in absolute measurements at the National Physical Laboratory [2], the primary winding 
consists of two equal single-laver coils of radius a, wound on the same evlindrical form se that 
each coil has an axial length of a2, with a yap between the nljacent coil ends of length a (see 
fig. 1, A) The two coils are jomed Ih series, magnet ally aiding The secondary ts a eireulat 
coil of square cross section, coaxial with the two parts of the primary winding and loeated in 
the median plane between their ends In this place, the racial components of the magneti 
field intensity, produced by the two primary windings separately, cancel for all values of racial 
distance, On the contrary, their axial components, which are equal, add, giving a resultant 
axial magnetic field intensity. With increasing radius this diminishes, becomes zero at: 
erttical radius, whieh is about 1.46 a, and for still larger values of the radius, ts in the opposite 
direction. The magnetic field intensity is, however, very small for moderate departures of 
radius from its critical value and for moderate displacements from the median plane. I 
therefore, the central filament of the secondary coil has a radius equal to the eritieal value, all 
the secondary turns lie in regions of small magnetic field intensity, so that the mutual inductance 
of the primary and a secondary turn is nearly the same for all the secondary turns, The 
resultant mutual inductance is closely N,.WM,. where MJ, is the mutual inductance of the primary 
and the central filament, and the small error may be accurately determined 

In the Wenner method for the absolute measurement of resistance [3] a modified form of the 
Campbell type of standard was employed. In this method the secondary is periodically short 
cireuited, making it desiable for this winding to have a minimum of inductance and resistance 


A coil arrangement was designed in which the number of secondary turns was only one-half the 





number on the primary This allowed the use of wire of relatively large cross section in the 
Present address, 10386 University Place, Schenectady &. N.Y 
Figures in brackets indleate the literature references at the end of this paper 
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Fiaurr 1 Vutual inductance standards 








A. Campbell type; B, National Bureau of Standards three-section type 


The necessary increase in the number of primary turns was obtained by the use of 
primary windings of greater axial length. As may be seen from figure 1, B, the primary 
winding consists of two end sections of equal length, and partly filling the space between then 
is a short middle section, The three sections of the primary are joined in series, aiding mag 
netically. The positions of zero magnetic field intensity were located experimentally by passing 
60-cyele current through the primary winding and using a short, suspended, magnetic needle 
whose period of torsional vibration was adjusted to the period of the current, to explore the 
field around the primary. With a chosen length of 5 em for the middle section and gap lengths 
of 5.6 em each, the field was explored for different lengths of the end sections. With relatively 
short end sertions two zero points were found with equal radii and located at equal distances on 
both sides of the median plane. With longer end sections these two points of zero field were 
replaced by two zeros in the median plane, locating two circular filaments of differing radii along 


secondary 


whose circumferences the magnetic field intensity is zero. The standard as used for the 
measurements had the following approximate dimensions: 

Pitch of primary winding 0.2 em Mean radius of primary 20.364 cm. 
Length of middle section 5.0 em. Mean radius of secondary 26.348 em 
Length of gaps 5.6 em. Number of primary turns 343 

Length of end sections 31.8 em. Number of secondary turns.. 218 
Approximate mutual inductance 10.897 mh. 


The radii of the circles of zero field intensity in the median plane were about 24.95 and 26.70 em 
The secondary turns were wound in a channel of rectangular cross section having dimensions of 
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2.542 em radial, 0.82 em axial, This arrangement of the standard will be designated as ease A 
in what follows 

For comparison, the dimensions of the Campbell type of standard, used by the National 
Physical Laboratory, follow 


Piteh of primary winding O.l em Radius of secondary 21.9 em 
Length of primary coils 750m Number of primary turns 150 
Length of gap between coils 15 oem Number of secondary turns INS 
Radius of primary Is oem 

Approximate mutual inductance 10 mh 

The secondary turns are wound in a channel of square cross section Lem ona side. It will be 


noted that the goal of a secondary of fewer turns, wound with larger wire, was reached in the 
NBS form 

In NBS Research Paper R P2029, page 296, the statement is made that consideration of the 
number and positions of possible points of zero field intensity is reserved for further study \s 
a result of correspondence with two of the authors of the paper on the absolute measurement of 
resistance by the Wenner method (James L. Thomas and Chester Peterson) the writer was 
encouraged by them to undertake a quantitative study of the NBS type of mutual-inductance 
standard, and they communicated curves of mutual inductance variations with secondary 
radius for two other arrangements with different lengths of end seetions. These will be denoted 
as cases Band C. In comparison with ease A, case B has two turns removed from the outer 
ends of each of the end sections, and case C has vet another turn removed from each of the outer 
ends of the end sections 

The manner in which the mutual inductance of a secondary turn in the median plane varies 


with the secondary radius is shown in figure 2 These curves are based on data supplied by the 
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National Bureau of Standards three-section type of standard 








authors of RP2029 and are reproduced with their permission.) For case A the mutual induc 


tance is a minimum for radius 24.95 em and a maximum for radius 26.70 em For case B the 
minimum is at 25.382 em and the maximum at 26.30 em. The curve for case C shows neither 


maximum or minimum in the median plane 


2. Method of Locating Points of Zero Field 
2.1. Zeros in the Median Plane 


For routine calculations of magnetic field intensity of solenoids, formulas previously published 
by the author [4] may be used; for definitive calculations the exact expressions in elliptic 
integrals must be used (see appendix for /7,). These formulas give the axial component of field 
/1,, produced by a current of | abampere flowing in a solenoid of radius a, winding density » 
turns per centimeter, and axial length 7 at a point at radius y in the end plane. Throughout 
this paper the symbol for current, 7, will be omitted in all formulas for field strength, it being 
understood that all values of 7/7 are per unit of current. Thus, for a point at radius y in the 
median plane, the axial field due to one of the end sections is obtained by making two caleula 
tions by the basic formula, one for length 2, and the other for length 2, (see fig. 1, B) and sub 
tracting the second from the first. The axial field intensity is the same for both end sections, 
and the values add in the resultant. The contribution due to the middle section is twice that 
obtained by making a calculation with the distance sz, (fig. 1, B Adding the resultants for the 
end sections and the middle section gives the resultant field, since, from svmmetry, the radial 
components of the end sections cancel, and the radial component due to the middle section is 
Zero 

The calculated values of //,/2n for the cases A, B, and C for points of different radii in the 
median plane are plotted in figure 3. Remembering the relation (d\//dy),~, — 2ry(/1,) 
it is seen that curve A checks the mutual-inductance curve figure 2, since the zero values of 
iH, oceur at y= 24.95 and y—26.70, the turning points of the mutual-inductance curve 
The maximum value of //, is found at y=<25.75, which is the point where the mutual 
inductance curve has its maximum positive slope. The curve B for case B (end sections 
shorter by two turns) is similar to that for case A; the zeros of 7/7, at y= 25.37 and y= 26.25 
agree in position with the minimum and maximum points of the mutual-inductance curve in 
figure 2, and the much smaller value of the maximum of /7/, is in line with the reduced slope and 
range of the mutual-inductance curve 

For case C, where three turns have been removed from each end section, the magnetic field 
intensity nowhere becomes zero in the median plane, although it is very small at y= 25.75 
This corresponds to the position of the minimum slope of the mutual-inductance curve in 
figure 2. There is no point of zero magnetic field intensity in the median plane for case C. It 
is, however, evident that, with a small fraction of a turn added to the end sections of case C 
the curve of /i, would become tangent to the zero axis at about y= 25.75, in which case there 
would be a single point of zero field intensity in the median plane at that radius 
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2.2. Zeros Not in the Median Plane 


The study of the mhagnetie field Intensity at pomts not mn the median plane Is « omplu ated by 
the necessity of calculating the radial component, //,, as well as the axial. Furthermore, the 
contribution of the end sections are not equal, and two distances « have to be considered for 
the middle section. ‘Thus, for calculating /7, and /7/, at a single point, the basic formulas have 
to be caleulated for SIX values of 7, both for I] and // For routine calculations, published 
formulas and tables may be used, but for final accurate values the elliptic integral formula for 
/7, given in the appendix has to be used. Then the radial component /7, may be made to depend 
on the formula for the calculation of the mutual inductance m, of two coaxial circular filaments 
one having the radius @ of the primary and the other that of the circle through the point 7? in 


question. “Two calculations have to be made for the two axial distances between /’ and the 


ends of the solenoid. If these distances are xz, and xy, with x, the smaller, then 
i 
// | m m, | 
2ry , 


For caleulating m, the Maxwell elliptic integral formula, or any suitable series expansion of 
this, may be used, or alternatively, the tables [7] published by the author may be employed 

The required zeros of magnetic field intensity are evidently the intersections of the loci of 
7, Oand 77,0.) The determination of the loci is made by calculating the field component 
for chosen points, for example, for a given value of the displacement & from the median plane 
and for a least three values of y. Thence the value of y for which the field component is zero 
is interpolated, From several interpolated points, so found, the curve of the locus may be 
drawn, When the approximate loci are not known beforehand, the exploration necessary can 
be quite time consuming. For greater accuracy in the location of a zero point, further caleula 
tions with the elliptie integral formulas may be made for the region of the graphically deter 
mined intersection of the loci, On account of the svmmetry of the loci about the median plane 
calculations have to be made for one half plane only 

In figure 4 are shown the loci and the zero points for cases A and C and for an arrangement D 
where the end sections are still shorter. The ratio of the length of an end section in comparison 
to the length of the middle section is 6.36 for case A; 6.28 for case Cy and 5.00 for case D 

The curves shown make clear the changes brought about by changing the length of the end 
sections with respect to the length of the middle section. In case D, the zeros are at radius 
24.7 and are displaced at equal distances = 2.7 em on either side of the median plane. In 
case C, the two zeros are at radius about 25.75, but are located only about 0.2 em either side 
of the median plane. In case A, the zeros are in the median plane with radi differing about 
1.75 em 

In all these cases the locus //,<0 is seen to have two branches, one in the median plane 
and the other, whose trace is practically the same curve in all these cases, intersects the former 
at about the same radius of 25.7 em. It is, therefore, easy to visualize the changes in the 
position of the locus //,—0 as the end sections are progressively lengthened in relation to the 
length of the middle section. With increase in this ratio, the intersections of the loct slide 
along the horizontel branch of the //,<0 locus until they merge in a single point at about 
y= 25.75 in the median plane, and, for still greater values of the ratio, two intersections again 
appear, now in the median plane, and these progressively separate, unsymmetrically with 
respect to y= 25.75, as the ratio is stil! further increased. The spacings are, for example, 

0.38 and +0.50 for case B; —0.8 and + 0.95 for case A. The proportioning of the lengths 
of the sections is quite critical in its effect. The striking difference between case B and C is 
brought about by the subtraction of a single turn from the outer ends of the end sections, 
with the gap lengths of 5.6 em, the length of middle section 5 em, and the radius of primary 
20.364 em, 


It is easy to show that decreasing the gap lengths, all else being unchanged, would shift 


the intersection of the two branches of the //,—0 locus nearer the surface of the primary, and 
thus with it, all the zeros to smaller radii 
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The negative slope of the mutual-inductance curves for the smaller values of the radius y 


in figure 2 raises the question whether there are still other zeros of field intensity, for the mutual 
inductance of a filament of radius y=0 is manifestly zero. The mutual inductance of a filament 


reaches in fact an absolute maximum for a radius equal to about the radius of the primary 
winding. In the median plane the component //, reverses sharply in passing through the 
winding, so there is here a zero of //,; but where no secondary can be usefully located 


3. Evaluation of the Correction for Secondary Distribution 


Even with the secondary coil placed with its central filament coinciding with the circle of 
zero magnetic field intensity, the mutual inductance will vary from turn to turn, For very 
precise measurements, it is necessary to take these variations into account 

If mo is the mutual inductance of the primary on the central turn of the secondary, and VN; 
is the number of turns on the secondary, the problem is to find the small corieetics factor 6 
in the equation 


M= Nym(1+8). (1) 
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The general formulas developed for 6 below vield a correction factor that takes into account, 
not only the distribution of the winding over the eross seetion, but also the effect <f any residual 
value of the magnetic field intensity at the center of the winding 


3.1. General Correction Formulas 


An obvious method of solution is to express the mutual inductance of the primary on any 
secondary turn by means of a Taylor's series development about the value for the central 
filament 

If my is the mutual inductance of the central filament, the mutual inductance m of any 
other filament, whose coordinates, referred to the center point of the cross section are (2), 
is given by 


dm dm | im d?m cm 
m— mote (Tr) aT) tafe (em) pee) +e (4) 


” 


| . d™m . d™m ' d™m 
ri eal Oo Raa] Ce Re 2 


the coefficients of the differential coefficients being those of the binomial theorem, The sub 
script zero indicates that the differential coefficients are to be evaluated at the central point, 

If the axial and radial dimensions of the rectangular cross section are, respectively, 6 and ¢ 
the total mutual inductance of the secondary is found by integrating (2) over the cross section, 
that ts 


V/ 7: midardy 


As the differential coefficients in (2) relate to the central point, they are constants, and the 
integrations are readily performed, and if it be noted that differentiation with respect to » is 
the same as differentiation with respect to y, the value of 6, carried to include fourth-order 
lifferential coeflicients, is evident by comparing the following equation with (1) 


‘ 1 (6? (dm = 1m bt ¢d‘m b*c? 7 d'm ce’ fd‘m ) 
M= Nim I ' mot24 ( dx ba dy ) ' 7900 \ da" ) y 576 (a dy) raa0 \ dy' )4 ia | 


Formulas for calculating the values of the differential coeflicients by two separate methods 
are given in the appendix. It should be noted, however, that certain fundamental relations 
exist between these coefficients 

If 77, and 7/7, denote, respectively, the radial and axial components of the magnetic field 
intensity, per unit current in the primary coils, then dm/dr 2ryll,, and dm/dy—2ry/1,, 
and in addition, the condition that curl /7—0 leads to the relation 


dH, dil, 
dr dy 


(3a) 


From these are readily found the general relations 

Pm l dm cm ) 
dy ydy dr 

d'm 1 dm d‘m_ . d‘m 

dy' y° dr’ dx* e dr'dy’ 

4 

d'm dm 20m, 30m 3dm f (4) 
drt dy y dy dy dy 
d*m 2dm 2m, 1m dm 


dedy yi dy yi dy ‘y dy dy J 
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correction 6 to be very 


(4) enables the second order terms by of the 
| 4 ° dm Co ‘bs 
Oe am | 24 ( dr’ ) ' o4y, | dy ) | 
l b?—e?* (d*m aes il \ 
“al oa Cae ) rae 5 | 


relatively very small, for the residual vaue of 
(Only when the 


Use of the first of eq 


simply expressed as 


The last term in (5) is the correction, usually 
the field 7/7, at the center of the cross section, due to imperfect centering 
mutual inductance of the central filament is a maximum or minimum will the value of be zero 

Use of equations one and two of (4) enables eq (3) to be expressed in terms of variable z alone 
There is found, however, no great simplification 


except for the terms in dm/dy and d*m/dz*dy 
National Physical 


In w paper on the calculation of a Campbell form of standard at the 
Dye [5] has made use of a formula for the evaluation of the correction for the 
Searle. No statement as to the manner of deriva 
It may be shown to be equivalent 


> Laboratory, 
distribution of the winding, due to G. F.C 
tion is given, and as printed it is not free from misprints 





Lo eg 3) as follows 


Substituting the first, third and fourth equations of (4) in (3) leads to the equation 


1 (1 /dm 4 b*¢ h' Pm 2 _ 5 4‘ h*¢ 
P ( ) "i ) : ) ' ( dy )C 24 ' 404, 2SSy ) ' 


: mol dy 24 mpg H40y5 
im he d‘m hb b*c c ) 
/ 
+ + 
\ dy )(. 57! 17) NOOY, 4 ( dy' )( 1020 576 aoe \ 
and noting that, in Dve's nomenclature, yo @ Ap, b= 2C, and ¢e 2B, and if we inelude the factor 
V, in each term, writing 7 in place of m to indicate this fact, we obtain Dve's expression free 


of misprints 


rr ” 20 B*C*—9C d 4 C? 9 C'—20 BC 
aM (G. | 6 360, } Cre) |= 6 * 360.4 } 


d in 10 °C? —6 C8 
Ci 360.4 ‘ 
To compute this correction, Dye calculated ddA for the whole secondary for a number of 


values of A, using the Campbell formula [1], which involves an elliptic integral of the third kind 
1, gave a nearly straight line with, of course, the value 


Peers 


These values of d\//dA, plotted against . 

of zero at the ordinate when the magnetic field intensity is zero 
Writing an equation for these calculated points in the form dM/dA 

e represents the difference between any radius A and that for which dM//dA is zero, the other 


51.6e+ 9.08, in which 


differential coefficients may readily be expressed, namely, d?M/dA* 516+ 18.06, and 
PM /d AS = 18.0 
For the center point, « 0.0210, so that (dMW/dA)y=1.138 wh/em. Other data were, 


Ap = 21.90, 2B=0.96, 201.00, and the calculated values of the terms in (6) are 
17.2 ppm of the total mutual inductance of 10 mh 
for which the 


a 15.008, 
0.00209 + 0.16894 + 0.00046 0.1715 wh, or, 
The value of the correction for this standard may be checked by formula (5), 


data become, a= 15.008, y= 21.90, ¢= 0.96, b= 1.00, Ny—485, and the winding density of the 


primary is #=—10 turns per centimeter 
From these there results (@’m/dz*)o=106.81 and (/7,)o=0.00204, and by formula (5 
0.024 16.946 ppm. The value of (/7,)) here found is less than 


the correction is b, 16.922 - 
This is explained by the 


the value 0.0169 derived from the value of dM/dA quoted by Dye 


fact that the radius of the secondary used was 21.8815 em, whereas the Dye curve showed that 


the radius for zero field should be slightly greater than 21.90, which was the value for which the 


calculations here are made 
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The fourth-order differential coeflicients were calculated by both methods A and B, (appendix 


which checked each other closely The values found are (d‘m/dr'), 5.381, (dtm /dr'dy? )o 
4.094, and (d*m/dy')o 2.584. Using these values in formula (3), the fourth-order correc 
tions terms are, in order, —0.1364-0.318—0.055— 40.127 ppm, and adding these to the main 


terms. the correction becomes 16.946+0.127 17.07 ppm, which ehecks Dve's value closely 
In Dve's formula, the first term; in formula (5), the last term, take into account the effect of 
the residual magnetic field at the center of the cross section It may be noted that, in each 
calculation, this amounts to only 12 parts in a thousand of the total correction for the cross 
section. The use of formula (3) with the equations in the appendix for caleulating the 
differential coefficients, involves materially less work than the Dye procedure 

The adequacy of the general formula (3) may be tested by calculating the correction for dis 
tribution of the turns of the NBS mutual inductance standard, case A. For this we have 
available the very accurately determined value, used in the experiments for the measurement 
of the ohm, NBS Research Paper 2029. This value was obtained from thirteen values of m, 
calculated by the absolute elliptic integral formula for equally spaced turns in the median 
plane. The effect of the smaller (axial) dimension of the cross section was evaluated expert 
mentally by directly measuring the difference of mutual inductance of two equal, flat coils, 
one placed in the median plane and the other at different axial displacements. The results of 
these tests, taken with the caleulated values in the median plane led to an estimation of the 


correction for distribution of turns as — 75.8 ppm, which may be considered as very accurate 


The data for the calculation by formula (3) are a= 20.364, yo 26.348, n-— 5, e = 2.546, 
b-O.82, 4, 2.500, ry, 8.100, and 2339.90. From the formulas in the appendix, we find 
I] ‘ O.046 (ad m dr ‘ 17 ov d'm dy' 16 0S, (dm dy dy 160] and ” pO.00] 


abhenrvs Whence 


l h? —« 2m . 
m ( 24 ) da ) S2.4 
| | 2xre : 
m, | 24 i ] t lo 
. aa de) | - OOK j P 76.3) ppm, 
Ny Q? dy 
| [ be dim ny 
m K da dy ) | =) + 4.5 J 
I e dém — 4 
m, EX dy' ) $7.0, | 


which agrees well with the assumed value, and indicates that terms of order higher than the 


80.8.) 








fourth are negligible 


3.2. Other Methods for Calculating the Correction 


It is evident that the process of computation is very time consuming when fourth-ordet 
differential coefficients have to be taken into account, because calculations of each differential 
coefficient, involving separate calculations for each of the three distances 7 have to be made 
It is, however, relatively simple to evaluate the second order terms in formula (5 The fourth 
order terms are, evidently, less important, the smaller the cross-sectional area. ‘This suggests 
assuming the cross section to be divided into two or more equal subsections and applying cor 
rection formula (5) to each. Thus the distribution correction can be found for each subsection, 


referred to the mutual inductance of its central filament. The average of the corrected values 
for the subsections, referred to the value of the mutual inductance of the central filament of 
the whole cross section, will give the desired correction 
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The details of this process may be illustrated for two subsections. For these we write 


\/ ~~ (Mo tf). \V/ 7 (My d 


\/ \ * z mm | \ mol | ead 


where my represents the mutual inductance of the central filament of whole seetion; mo, and 


and, therefore, 


Mow, the values for the central filaments of the subsections, and the total correction is 


A Miay ot My A, ] d d 
sn +— , » with 3 [” a : m ] mo and a ( = ) 
é Mig 


If the subsectional area is sufficiently small, only the second order terms formula (5) need be 
calculated, but the three mutual inductances my), mo, and my will have to be calculated by the 
precision elliptic integral formula 

(s an example, assume the NBS secondary, case A, to be divided into two subsections, 
centered in the median plane. For each subsection ¢= 1.273 and 60.82. Taking the radii 
of the centers of subsections as y,, and y, and my—49991, as before, the details of the caleu 


lation are as follows 


- Mio Mhy & 
y 25.711 Y. 26.985 a 199OS8S8.05 
d mm a Pm oe ous ><> 
( - ) L276 ( dr? ) 28.288 m, 19001 32 
// 0.0704 (77,), 0.03608 Difference 3.27 
d 0.0849 d. 1.1376 A 65.5 ppm 
d,+d 
-s 0.5263 Meo) 19084.20 a: ros 10.5 
4 J a 10.5 ppm Mos 190901 .9] 6 76.0 ppm 


Evidently, in case only the correeted mutual inductance is required, it is only necessary to 
calculate M/, and 7, and to take their sum. This suggests still another procedure, using the 
sectioning principle 

The secondary may be replaced by two equivalent filaments, according to the method of 
Lyle [6], and the required mutual inductance, corrected for the distribution of the winding, is 
approximately the average of the precise values of the mutual inductances calculated for 
these equivalent filaments 

The radii of the equivalent filaments are given by (r—D) and (r+ D), where 


} wt og (Om) and J)” 


12 


For case A, r= 26.54906 and ))= 0.695804, so that the radii of the equivalent filaments and the 
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mutual inductances, corresponding are 


y, =r—D=25.65326 m 1998336 
Yo —=r+ D=27.04487 m 19901 .42 


Average — 49987 .39 


m LOOT 2 
Difference 3.98 abhenrvs 


78.6 ppm 


Lyle’s formulas neglect differential coeflicients higher than the second, so as would be ex 
pected, this result for the undivided secondary refleets the neglect of fourth-order terms low 
ever, the secondary may now be assumed divided into two subsections as in the previous 
example and the two Lyle filaments for each of the subsections considered 

This treatment yields the following results, the values Am, referred to mo. being given instead 


of actual values of m: 


Y= 2Z5ABISI Am, 10.5387 
Yo 2599368 Am SASLS 
y °26.70446 Am L505] 
Yq 27.26668 Am, 2.5946 
Average 3.8024 abhenrvs 


The distribution correction resulting — 76.06 ppm 

Evidently, by continued subdivision, with separate treatment of each subsection and com 
bination of results, any desired accuracy may be attained. The labor, however, increases 
rapidly with the number of sections. The check furnished by the known value of the corres 
tion for the NBS coil, case A, indicates that the caleulation using the second-order terms, 
using two subsections, or the calculation of the average of the mutual inductances of the four 
corresponding Lyle filaments is sufficient for practical purposes, even in the most precise 
work 


3.3. Examples of Special Cases 


From formula (5) it is evident that for a square cross section the first term is zero, and if the 
central filament is accurately placed at the position of zero magnetic field intensity, the rela- 


tively small second term of (5) is zero also. That is, the distribution correction in such a 


case depends entirely on fourth- and higher-order terms 

In the case of the standard calculated by Dye, the correction depended largely upon the 
factor (B°—C"*), which had a small value due to a small departure of the cross section from 
being a square. However, assuming an accurately square cross section, accurately centered 
on the radius of zero magnetic field intensity, the correction for the Campbell type of standard 
should be very small. Also, that arrangement of the three-section standard, denoted by the 
name, case C, may with proper design be made to depend upon the theory for a square cross 
section 

For a more detailed study two cases were selected: (a) a Campbell standard with square cross 
section secondary; and (b), the arrangement in figure 4, case C. In order to obtain comparable 
results, cross sections of nearly equal areas were assumed in each case. Case A, already studied, 
has a cross-sectional area of 2.09 cm*. For case C, a section with an axial dimension of 2 em 
and a radial dimension of 1 cm, area 2 cm’, was assumed, and for the Campbell type a section 1.5 
em square, giving an area of 2.25 cm* was assumed. 

A detailed study of the Campbell case is covered at length in section 4. The correction for 


cross section is found to be about 0.70 ppm. A separate calculation by the general formula (3) 
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vave correction for 


Residual field O.OLS 
fourth-order termes 0.282 
Total 0.245 ppm 
For case CC) using « 1. d=2. » 95.75. (H 0.00304, and (d?m da O.7153 there os 
found 
Second-order terms LSI 
fourth-order terms 1.15 
Total 0.00 ppm 


On aecount of the symmetry of this case with respect to the median plane, a more accurats 
method of treatment is to consider either square half of the cross section. For these the residual 
field intensity at the center is (H 0.0234, so that we have by (5) the correction 0.123. The 
fourth-order terms amount to 0.323, so that, referred to the center of the square the corres 

tion is — 0.20 ppm But it is found that the mutual inductance of the filament at the center of a 
quare is —0.42 ppm referred to the center of the whole section. Therefore, referred to this 
latter point, the correction ts — 0.20 plus — 0.42, giving the result 0.62. which agrees well with 


the previous determination 


4. Detailed Study of Variations of Mutual Inductance From Turn to Turn 
Over the Secondary Cross Section 


The variations in mutual inductance from turn to turn of the secondary coil can be obtained 
by direct calculation of the mutual inductance for filaments with regularly spaced values of the 
radius y and the displacement & from the median plane As the value at each point ts the 
resultant of caleulations with six different values of 7 in the fundamental elliptic integral formula 

for the NBS type) this is very time consuming, and on account of the near cancellation of 
terms, the resultant is appreciably affected by small errors in the individual terms 

It is more logical to make caleulations of the sought-for changes AM directly, than to obtain 
these by subtraction of calculated values of 7, and greater accuracy in the result ts attainable 
For this purpose, two series expressions for M are available, one suitable for larger values of 4 
due to Rosa [8], the other, suitable for smaller values of 2, given by Dwight [9]. Noticing that 
for a point 7? having a displacement from the median plane, referred to a point having the same 
value of y, but situated in the median plane, a calculation of twice the value, using the longet 
distance 2,, has to be made for the latter, whereas, for the point ?, the contribution of the left 
end section is made with the distance (2, + &), and for the right-hand end section with the distance 

r,—é Combining these, a difference formula may be derived for the change A’.\7, due to 
displacement, as far as distance 2, is concerned. Similarly, the application of a difference 
formula for distances 2) and zy, enables A’’M and A’’’M to be found. The desired total change 
isgiven by MAW AM AU MY AM. This change A.V gives the difference in mutual induet 
anee of point /? from that at the point in the median plane having the same radius 

Likewise, these two basic series formulas may be applied to find the difference in mutual 
inductance between a chosen point in the median plane, having a radius y, and the reference 
central filament of radius y%. These latter calculations enable us to refer the calculated AA/ for 
all points to the reference value of mutual inductance of the central filament. 

Details of these difference formulas are given in the appendix. Even using the difference 
equations, the results have still to be found by combining terms that partly cancel, and careful 
checking is necessary to detect errors. The work is necessarily arduous 
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5. Distribution of Flux Over the Cross Section 


5.1. Campbell Type 


For purposes of comparison, a primary having the same winding density and about the 


sume radius as the National Bureau of Standards type is assumed, the details of its design being 


Radius of primary —a— 20 em 

Winding density, n= 5 turns per centimeter 

Axial length of each primary section — 10 em 

Mean radius of the secondary — 1.4582 a— 29.164 em 


With this secondary radius, the secondary turn at the center of the cross section lies very closely 
in the circle of magnetic field intensity zero 

There are LOO primary turns and 594 secondary turns: the mutual inductance would be about 
the same as that of the National Bureau of Standards mutual-inductance standard 


Assuming for the secondary coil a square cross section, 1.5 em on a side, 49 values of mutual 














nductanee value referred to the value at the center of the cross section. were ealeulated 
wing the difference formulas of the appendix The grid of these value Which are spaced 
at intervals of 0.25 em of Ay and 0.25 em of distance & from the median plane are shown in 
firure 3 The maximum deviation from the value at the center nbout S50 ppm Phe value 
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eau arée.5 710.41 ree6e ro rai oh 
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L l l | l rt 
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t 
Fieurr 5. Calculated values of variations of mutual ind ” 
rom turn over the cross section Clam ple fype 


of MMV are symmetrically arranged with respect to the median plane and are mostly plus 
those in the median plane are nearly symmetrical about the mean radius and are negative 
It wes noted, during the calculations that, for a given radius, AM is very closely proportional 
to & This is true also for the other cases investigated 

Wilmotte [10] has derived expressions for the change 6.7 in mutual inductance for point 
of the cross section of the Campbell type of standard, referred to the value at the central 
point of the section His eq (4) expresses this as a function of the coordinates x and y, referred 


to the center point as origin, and differential coefficients of the radial component of the mag 
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netic field. This expression, which has recently been brought to the attention of the author 


in form seems to be similar to the formula (5) of this paper 

Wilmotte’s eq (5) expresses 6.17 in polar coordinates (47,0), referred to the center of the eross 
section as origin, The principal term of this shows that 6.17 is proportional to dr cos 26 An 
inspection of figure 5 shows that this law is here approximately obeyed. The slight lack of 
ymmetry about the mean radius in figure 5 is in line with the second term of Wilmotte’s 
formula (5), which is proportional to sin’@, and involves higher-order differential coefficients 
of the radial field component. This term would add to the principal term amounts of opposite 
“yn for points of greater radii than the radius of the center point than it would for points of 
smaller radius than this value. Tt may therefore be concluded that formally, at least, th 
results of this investigation confirm Wilmotte’s equation 

To obtain the average deviation of the mutual inductance of a secondary turn from that 
of the central turn, the deviations have to be integrated over the cross section and the result 
divided by the area of the eross section, This integration may be accomplished by means 
of a suitable mechanical integration formula, such as Simpson's rule Using this with 
the values of AA of figure 5, there is found a correction for cross section of + 0.71 ppm 

A check on this value is provided by the well-known quadrature formula derived from 





the Taylor's series formula (4). This gives the average mutual inductance of an array of 
filaments distributed over a reetangular section 2Ay by 2& in the form f 
; 
\ 2 
\V/ (A0,+ M+ M+ A+ 2M), 7 ) 
t 


in whieh 


Vi, the mutual inductance of the central filament, 
VJ, the mutual inductance of the filament at (Ay.0 


SET pe 


V/,— the mutual inductance of the filament at (—Ay,0), 
VJ, the mutual inductance of the filament at (0. 
V/,- the mutual inductance of the filament at (0,—&) 


Replacing these by My w+ d&. Wy) wt}, ete., in which » is the mutual inductance of the 


reference filament, which ean be outside the area considered, we have the formula 


VW Nowl t+ § CAAT + Ay AD Ag AD AGATE 24,41), (8 


gE + 


which expresses the average mutual inductance for any subdivision of the secondary cross 
section in terms of the mutual inductance » of the central filament of the secondary and the 
deviations AV of the grid of values. Applying this quadrature formula to the nine component 
square sections of which the secondary may be considered to consist, and remembering that 
each seetion consists of only one-ninth of the total secondary turns N,, we have to average 
the following nine values 
25. 18, 855.18, Y51S, 
365.08, —0.13, +365.08, 
} YS.04, 374.38, 4+ 28.04, 


which vields the total correction for cross section of + 0.69 ppm, which checks the other inte 
gration method, The correction was calculated also from the Taylor's series formula, with the 


result 
Correction for residual magnetic field at the center +O, 018 
Fourth-order terms +0, 232 
Total ppm +0, 25 


The reason for this surprisingly small correction may be made clearer in what follows 
In the median plane the mutual inductance is a maximum for the central filament. For 
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ay 


a given radius, the mutual inductance is a minimum in the median plane. Thus, 


the pout 
of zero magnetic field intensity is located in the center of a saddle-like region of flux distribution 


This is made clear by figure 6, which is a plot of the loci of points having the same values of 
VM. These curves are obtained from graphical interpolation from the grid of caleulated 
values of AM in figure 5. In figure 6 each curve is labelled with the corresponding value 


of AMZ in ppm. Each curve is the locus of filaments linked with the same flux. That is. each 
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locus traces a line of magnetic flux, and the figure gives a pieture of the flux distribution over 
the eross section of the secondary Assuming the positive direction of a flux tine along the 
axis of the primary as from left to right, the direction of the lines in the lower quarter of figure 6 
is from left to right, and of those in the upper quarter from right to left. Those in the right 
hand quarter are directed inward toward the axis, and those in the left-hand quarter outward 
from the axis 

Figure 6 makes clear that the very small total resulting correction for cross section is the 
result, not of small deviations of mutual inductance of the filaments, but from the very close 
balance of positive and negative deviations 


5.2. National Bureau of Standards Type, Case C 


The treatment of this case follows the same lines as for the Campbell type. A grid of 45 
calculated values of AAS spaced at 0.25-cm intervals of Ay and 


by Lem in dimensions, is shown in figure 7, These values of A.M are distributed symmetrically 


£ overt the cross section, 2 em 


about the median plane, as would be expected, but about the radius of the center of the cross 
section (25.75 em) positive values on one side are nearly offset by negative values on the other 
Furthermore, excepting fer the corner points, the actual magnitudes of the AM are less than 
100 ppm, compared with values 5 or 6 times as great in the Campbell type 


For a given radius, AM, is closely proportional to &, a fact that does not appear in the 


figure, where all values are referred to the value for the point at the center of the cross section 
In figure 7 are sketched in the loci 17,<0 and 7/,—0. Along the former dM/dt—0 
and along the latter dM//dy—0. These facts are checked in a general way in the figure 
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Averaging the values of figure 7 over the cross section by Simpson's rule, the correction 
found is —0.44 ppm, and by applying the quadrature formula to the four square sections 
into which either half of the total cross section may be divided, the result is —0.38 ppm 
Either half of the cross section may be replaced by the Lyle filament to which the square 
section is equivalent. This has a radius of 25.75[1 4 (1/24)(1/25.75)*], and the coordinate 
‘ 


x - 


0.5. It is therefore, closely coincident with the filament, (25.75, 0.5) in figure 7, for which 
AM 0.42. Thus this value may be taken as the value of the correction given by the Lyle 
method. Thus, the conclusion is amply checked that, for case C, the correction for cross 
section is not far from —0.4 ppm 

Figure 8, which is derived from the values of figure 7, gives a picture of the distribution of 
AM over the cross section and the shape of the flux lines over the whole region. This figure 
demonstrates in a striking way the smallness of the deviations in mutual inductance of the 
filaments of the secondary. The directions along the flux lines are indicated by arrows on 
the loci of equal values of A.V, and also there are indieated the loci of //,<0 and 7,0. At 
their intersections with the curves for //,<0, the AM should have a radial direction, and this 
is substantiated by the curves as well as can be expected. (C, and ©, are the points of zero 
magnetic field intensity 

It may be noted that, in passing in a circuit around the center of the cross section, there 
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are six regions of the field in which M47 is alternately plus and minus, whereas for the Campbell 
type there are four This suggests that with a primary of » sections there would be 2n such 
regions 


5.3. National Bureau of Standards Type Case A 


For comparison with ease C, figure 9 has been plotted for case A. The data for this was 


obtained for the median plane from the curve of mutual inductance for points in the median 


> 


plane, figure curve A. Values of the difference in mutual inductance at the boundaries 
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3 0.41 and for the median plane, at the same radius, were found for three values of the 


radius, namely, in the center and at the extreme corners. These values are relatively small, 
ranging from +59 ppm for the larger radius, + 29 for the mean radius, and — 47 for the smallest 
radius. From these, other values on the boundary were derived by graphical interpolation 
For interior points, the values were obtained from the law of proportionality to &, which 
applies in all these cases (see fig. 19, RP2029). The grid so obtained checked the axial cor 
rection of + 7.2 ppm given on p. 316, RP2029, to one-tenth part in a million 

In figure 9 are sketched also the loci //,—0 and //,—0. The points of zero field intensity 
are indicated by ©, and ¢ The directions of the flux at different parts of the field are indicated 
by arrows on the lines of flux mapped by the loci of equal values of AMZ, and also in the margin 
of the figure are indicated the directions of the axial and radial components of the field) in 
different regions The corresponding diagram to the left of the cross section are not indicated 
but they are at once evident if itis remembered that the radial component reverses in passing 
the median plane, and the axial component is continuous 

The range of variation of mutual inductance over the cross section in only about 200 ppm 
on either side of the average, which is nearly as good as for case C, and 4 or 5 times better than 
the Campbell type, but the marked lack of symmetry about the radius of the center of the 
cross section is evident. The center of gravity of the values (the point where \/ 75.8 
ppm) lies at about y— 25.97 em in the median plane, that is, about 0.55 em below the center 


of tie cross section 


6. Conclusions 


1. The inherent balancing out of the plus and minus variations of the mutual inductanes 
from turn to turn over a square cross section has been noted and illustrated. This is in line 
with the well-known fact that, in mutual-inductance calculations, the geometric mean distances 
of a square area from an external point is very closely equal to the distance from the center 
of the square to the external point 

2. In the Campbell type of standard, for this reason, the correction for cross section is 
surprisingly small, even with rather large variations of the mutual inductance among the turns 
ver the cross section, Llowever, in order to be assured of this satisfactory balance, uncer 
tainties in the actual placing of the turns must be minimized. The example here illustrated 
assumed a cross section greater than is necessary. The NPL standard employs a square 
cross section of secondary only | em on a side, and, for such a coil, the extreme variations 
AM would be of the order of less than half those in the example above. Where the use of a 
finer wire on the secondary offers no disadvantage, the Campbell type has the advantage of 
simplicity 

}. For a given area of secondary cross section, the three-section standard has been proved to 
have much smaller variations 4.1/ over the cross section than the Campbell type. This is true 
both for case A and case C. However, case A shows a lack of symmetry of the distribution 
about the central point, and the balancing of values of opposite signs is not taken advantage of 
The correction for cross section may be calculated to better than 1 part in a million, but its 
actual value appears large in comparison with the other cases 

An adjustment of the shape ratios of the three-section standard so as to attain the special 
case of a single circle of zero magnetic field intensity in the median plane, and the use of a second 
ary of square cross section, accurately centered in this zero region, would seem to be an optimum 
arrangement, but the adjustment would be very critical. A good approximation to this is 
offered by the arrangement case C. In the example of this given above, there is a secondary 
that is essentially a combination of two equal square cross sections of equal mean radius, and 
the symmetry about the median plane leads to a total correction for cross section which is the 
same as that for either half of the cross section. Furthermore, each square component is cen- 
tered in a region of nearly zero magnetic field intensity. The range of the values of A. is in 
general much smaller than the extreme values of about +150 ppm, and it is probable that the 


secondary could actually be wound to have a correction for cross section not exceeding a part 
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ina million \ refinement on this case would be to slightly adjust the dimensions so that the 
two zero points would be located exactly at the centers of the component squares that form the 
cross section, but this would also be a critical adjustment 

1. A thorough study of the distribution of the AA over the cross section by ealeulating the 
values for a grid of points is very time consuming. The abbreviated methods already discussed 


will. in general, be used. The Tavlor’s series formula should be applied to subsections of di 


mensions sufficiently small to require the caleulation of second-order terms only. However 
if the field intensity //7, is not zero at the center of the subsection. its value must be caleulated 
and correction made. For case A, two subsections are found to sufliee. Since. however 
caleulations of (// and (d?m/da are necessary for each, the caleulation of AZ directly for 


equivalent Lyle filaments is to be recommended. Especially is this true, if the given cross 
section may be divided into square subsections. For each of these a single Lyle filament is 
required and the average of the mutual inductance calculated for each of these, averaged to 
find the final corrected value of the mutual inductance. It is only necessary to be sure that 
subsections are of small enough area to ensure that fourth-order differentials are negligible 

With case A it has been shown that division into two subsections, not squares, calculating 
and averaging the mutual inductances of the four corresponding Lyle filaments gives a very 
accurate result It is interesting in this connection to note that the cross section in this case 
may be divided quite closely into three square subsections. If this were accurately the 
cause only three Lyle filaments would be hecessary, and the subsections would be so small as 
to make certain that fourth-order differentials would be negligible 

For case C, we have an especially favorable case with the cross-sectional dimensions chosen 
i. @., O=2, « 1, in that the Lyle filaments for the whole cross section are the same as the Lyle 
filaments for the two equal square subsections. The residual field is verv small over each sub 
section, and higher-order differentials are negligible. The caleulation of the mutual inductance 
of the single Lyle filament for a square subsection gives a very accurate determination of the 
corrected mutual inductance per turn. The final value to be used for the mutual inductance 
of the whole seecundary is merely NV, times this value 


7. Appendix 
7.1. Evaluation of Differential Coefficients 


Methods for calculating the differential coefficients that appear in the general Taylor's series 
formula for the evaluation of the correction for the distribution of the winding may be based 
on two separate methods of approach, (a) on the formula for the radial cemponent, per unit 
current, of the magnetic field intensity of a circular filament exerted at a point 7? distant y 
from the axis and » from the plane of the filament, and (b) on the formula for the mutual 
inductance of two coaxial circular filaments, m,, one of them of radius a and the other y, the 
distance between their planes equal tov. Since in each case computations have to be made 
for the three distances x, 7, and vy, the basic elliptic integral formula has to be used in both 


methods, no single series development being sufliciently convergent for all three values of 4 


a. Method A 
Placing //, for the radial component of the field, due to a unit current, ina cireular filament 


of radius a at a point 7’ distance y from the axis and s from the plane of the filament 


(2) — —2ey(2n) es —Hor+ Hes . 


Calculations made for the three distances, 7), 2, and z, are combined as shown; » equals the 
winding density of the primary winding, and the factor 2 takes account of the symmetry 


about the median plane. 
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1, is caleulated by the Russell elliptic integral formula [11] 


Ho [2 BKB) } 
ry ’ 


ire the complete elliptic integrals of the first and second kinds to modulus 4 


in which K and Ie 


The coefficient (d?m dy 


is the residual 


in which (// 


center of the cross section 


In this expression 


Numerical values of A 


of reference [3], or alternatively, 


e’kr| ¥, where ¥ 


and /, 


lely a) 


of vy, A, may be 


axial component, pet 


IT is an ellipti 
I} may be computed by the methods illustrated in tables 5, 6, and 7 


y is Legendre’s elliptic funetion defined in eq (28 
calculated by 
\5) 7/7, must be evaluated for the three distances 7, combined as in eq 


is derived from the relation 
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For the fourth-order differential coefficients we have 
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in (A6) has to be calculated for each of the distances 2 


The common factor 2n takes into account the winding 


of the primary and the symmetry about the median plane 
The differential coefficients of //, were obtained by differentiation of the elliptic integral 


formula (A2) 


of ways in each case 


where 





dH 22 wl K—E 2ay ok | 
dy ry ? 


y . 2ay*(y—a 
rst a, 


a=) 
r rr3 
y(y+a) , 2y(y—a) 
} meee t : 
r r 
CH, 2x . , 2a : 
<= | (K—E)—“"” cE |, 
Yd ny r 


316 


A2 


A3 


\4 


Lat the 
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This offers no especial difficulty, but the result may be transformed in a number 
The forms which follow are believed to have some advantages 
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The modulus of the elliptic integrals and the equations for 7? and 72 are the same as in the 


i formulas (A2) and (.A5) above 


b. Method B 


‘mpen 


In this method of approach, the base is the formula for m,, the mutual inductance of two 


oaxial circular filaments of radi a@ and y, with a distance x between then planes 
: dem dm dm dm : 
H ( da ) nan ( da ) ( da ' ( da | AIS 


which indicates that calculations have to be made for the three distances w,, a, and vy and the 


aN pe Sree ny 


\ results combined as in (ALS The common factor 22 takes into account the winding density 
of the primary winding and the symmetry about the median plane 


The fourth order differential coefficients are 


iq dtm a dim im im 
. ( dr ) - ( ds ( da ) +( da ) | 
\l4 
i dim " im dm dim 
(dy ) 7” i (Su) ( inde) He) | 
; and from these and the general relation 
dtm | dm d*m e dim J 
( dy' ) y ( da ) da ~ dardy als 


is derived the value of (d*m/dy' 
Starting with the Maxwell elliptic integral formula for the mutual inductance of two coaxial 


circular filaments in its integral form 
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the required differentiations are made under the integral sign. In the results occur integrals 
of the form 


"2 de 
/ | (1—k* sin? gg)" 
which are related by a recursion formula. We have, in fact, 
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Introducing then the nomenclature 
pla Pl,—@l, pep ls—_€l,, AIS 
in which p 9— kk) and @-2/k*: the differential coefficients may then be written in the 
forms convenient for calculation 
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In all these elliptic integral equations the modulus and nomenclature is the same as above 
For each differential coefficient calculations for all three values of « have to be made and com 


bined in the standard way 


Results of caleulations of numerical problems give results by the two methods, show that 


methods \ and B closely ayree in all CUSCS 
7.2. Difference Formulas 


These give the difference in mutual inductance AA between the value for a circular filament 


in the median plane and two equal solenoids of length 2, and the value when the filament is 


displaced a distance & toward one solenoid and the same distance away from the other 


For the longer distances, «,, Rosa’s formula [S] is made the base 
V/ Jr aapnsr@ 
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+q, and (+ dq are derived from the expression for Q by replacing x, by 
respectively, and p by (p+r,) and (p47 
Finally, 4,17 ts obtained from 3,.47/\7, using equation (a) to caleulate MV 


For shorter distances u ora the basic formula adapted from Dwiwht 1] Is 


Vi 2a ryauG 


in which 


with 


The distances u Is, 35 tf} are substituted for x to find G,@ 4), and @ 


AM 
V/ 


from whenee, using r=, in (b) to caleulate MV, we find A, 
To caleulate the difference between the mutual inductance for a filament with radius y, 
situated in the median plane and the mutual inductance of the central point with radius y, 


we have for the longer distances 4 


AM ; a (O&O ss q | 
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in which & yy) Gf +45) and ( is calculated from the formula for Q@, with 2 .2,, and with 
y= Yo in the equation for p, and (4+ q) is the corresponding value from the formula for ( 
with rv, and y used in the equation for p 

For the shorter distances xz, (or 2), the formula becomes 


AM € € 5 q 
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in which €= (yy) Yo and Gy is derived from the above formula for G with #— 4, and yo my, 


and g is the increment in this when y is replaced with y 
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‘The Deuterium-Sulfide Band at 4,590 cm 


Harry C. Allen, Jr.,* Robert E. Naylor,’ and Earle K. Plyler 


The absorption of deuterium sulfide has been observed and measured under high 
resolution in the region from 4,513 em=! to 4,675 em rhe rotational structure has been 
analyzed through the use of published energy and line-strength tables for the rigid asymmetric 
rotor A classical centrifugal-distortion correction was applied to the rigid energy levels 
Che absorption is attributed to the vibrational band having an excited state with the quantum 

umber (m;,%,n 11,1 The inertial parameters giving the st fit to the observed 
absorption are 
Ground stats 15.473 em Bh 4.490 2.444 
Excited stats l 5.465 en B 1.482 ( 2.383 
L sing the ground state inertial constants, the D-S distance is calculated to be 1.345 A. and 
the D-S-D-angle 02°16’, in excellent agreement with similar structural parameters determined 
for hydrogen sulfide from infrared studies 


1. Introduction multiplier, which detected the fringe system, and 
the other circuit was connected to a PbS cell suit 


able amplifiers were used with each detector. Stand 
ard atomic lines from krypton were thrown on the 
fringe system, and in this manner all the maxima of 
the fringe system could be reduced to wave numbers 
The instrument was run continuously while the 
standards, fringe system, and absorption spectra were 


Although the rotational fine structure of several 

red bands of hydrogen sulfide has been satis- 
factorily analyzed {1 to 6], there has been no report 
of a successful analysis of the rotational fine structure 
of any infrared band of deuterium § sulfide. The 
spectrum of DLS has been observed under rather low 
esolution [7, S], but no satisfactory analysis of the 
otational structure was presented In order to 
nable a solution for the cubie and quartic terms of 
he vibrational potential function of H,S, the infrared 
pectrum of DLS has been reinvestigated under high 
esolution 


recorded 


3. Vibrational Assignment 


The absorption at 2.2 « shows well-defined 7? and 
PR branches with an extremely strong absorption near 
the center, which suggests a collected @Q branch, 
This absorption corresponds to an energy range for 
which there are several possible vibrational assign- 
ments. Dennison [10] has shown that as p=B A 
becomes smaller, the Q branch of A-type bands tends 
2. Experimental Procedure to col'set at the center. On the other hand, there 
is a gap at the center of a B-type band for all values 
of p. Henee, it is concluded that this band is an 
\-type band, i. e., the electric moment oscillates 
along the least inertial axis. This axis is perpendicu 
lar to the symmetry axis; thus, m,; must be odd, leav- 
ing only two possible vibrational assignments (7 
No, 1 11,1) and (0,3,1 Arguments relating to 
the relative intensities, to the dependence of the 
spectrum was obtained with a pressure of 42.5 em | Moments of inertia on n2, and especially to the value 
f He. In this wav the weak and the strong lines of A=1-—T,—T1, indicate a strong preference for the 
ould both be measured to the best advantage. The | fetmer assignment. However, the rotational analy- 
spectrometer had a grating with 15,000 lines per | 5! is independent of the vibrational assignment ex- 


uch, and the resolution was about 0.10 em-'in the | CePt for the effect of the symmetry of the excited 
1 revion vibrational state on the rotational selection rules 


In this paper an analysis of the absorption of D,S 
the region from 4,513 em™! to 4,675 em™! will be 
press nted 


The intensity of the absorption bands of DS in 
i region from 1 to 3 «4 is low. With a cell of 6-m 
path length it was possible to observe the spectrum 
of the 1.1.1) band of D.S in the region of 2.2 u 
vith a pressure of a few centimeters of Hg. The 
spectrum of this band was measured for several 
pressures, up to 60 em of Hg. In figure 1 the 


In order to measure the lines to a high precision a 
Fabry-Perot interferometer was employed in a man- 
her similar to that used by Douglas and Sharma [9] 
\ two pen recorder was used for observing the spec- 
trum. One cireuit was connected to a 1P28 photo- 


4. Analysis 


Initial trial values for the inertial parameters and 
the band center were estimated from the known 
parameters of hydrogen sulfide [11]. A trial spee- 
work wes eganeneed bey the Atemte Bauseny Conmigo trum was caleulated through J=8, using published 
t address: Michigan State College, East Lansing, Mich tables of asvmmetric-rotor energies {12}, and using 


vidress: Harvard University, Cambridge, Mass 
brackets Indicate the literature references at the end of this pa per A-ty pe selection rules 1. @., dipole change along the 
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A classical centrifugal-distortion 
rigid energies The 
intensities were calculated with the aid of 
published tables of line strengths [13], using the line 


lead inertial axis 
correction [1] was applied to the 
relative 


strengths tabulated for K=0.5. These line strengths 
were combined with the Boltzmann factor and the 
appropriate nuclear-spin factor For DS the syvm- 
metric rotational states are twice as numerous as the 
antisymmetric rotational states, hence the relative 
intensities of all transitions in which the ground state 
ris even are multiplied by 2. 

The calculated spectrum was compared to the ob- 
served absorption. The adjustment of the six iner- 
tial parameters and the band center to give the best 
fit to the observed spectrum is largely a stochastic 
procedure. The inertial parameters were adjusted 
by means of the partial derivatives of the energy 
with respect to the initial parameters [1] until it was 
possible to assign unambiguously 100 transitions for 
J<S8. In some instances two or more transitions 
were assigned to the same peak. Final values of the 
inertial parameters were determined by solving 100 
equations of the type 


E wn EF eax An +a’ of +h’ Of 
OA oR 
(7 QE _ dE ,dK_ dE 


se U 7: 
oC OA oR “or” 


by the method of least squares. The best values of 
the seven variables are given in table A compari- 
son between the observed absorption and the spec- 
trum calculated from the constants given in table 
is shown in figure |. In the final calculations the 
mean deviation of about 110 transitions, with J< 10, 
was 0.08 em*'!, 

An inspection of figure | reveals that most of the 
absorption in the central region is well accounted 





The correction may be readily evaluated for individual energy levels, and 
hence has this advantage over methods that involve tedious calculations 


‘ fy 
Aim nahi i A 
Wu WwW | WiiMA MMW hile \ 
vey ‘V 
‘ 


a ~ 
Ar Lorhbord Uw 
ve Vv wy WW VV 





' 

? F > PP de obsess vee aoe « 
se ciddddddcdisddddad ds g 
eee ee *> > eeweeoe = oh ae wt ' 
SCS PFAeGerrdéddddddnanan IF a 


Fieurs 1. 


AWW 


P 


eS Cee vs i. a re wii aid Juhl 


yeas 28 Nand dodsdr ow 
Toor 75 ? awk ake, 
nate ~ a n'e nny enh worm aed 
ten aw Ad dddodrnd 


The infrared-absorption spectrum of deuterium sulfide in the region from 4,500 to 4,700 em 


Tarte | Constants of the DoS molecule 
r EF xvcited 
i 473 em 65 en 
B 141m 14s 
( its ist 
/ 114 121 l 
j ¢ a) 6. 244 
/ 11 451 11744 
Bond ang) "le 
D-S distu 145A 
1,11 my $2 on 
for. Some of the weaker peaks are due to a piling 
up of weak transitions that were not plotted. A 


the sample of D.S contained about 5 percent of 
HDs, and some of the weaker peaks may arise from 
absorption by this molecule. The absorption or 
the ends of the region that is not accounted for arises 
from higher J transitions. Transitions for J > \2 
were not calculated because they are beyond thy 
range of presently available energy tables, and the 
labor involved makes the task unrewarding 

That some errors between calculated and observed 
values are greater than the mean deviation jis 
ascribed, mainly, to the effect of a finite slit on 
several close-lvying transitions, the integrated result 
not corresponding uniquely to any single transition 
Likewise, some error is introduced by the use of a 
classical centrifugal-distortion correction, although 
it is believed that this will not cause errors greate: 
than 0.20 em™' in the most unfavorable 
Table 2 gives the assignments and calculated lin 
positions for the main P and F& branch transitions 

For a nonlinear triatomic molecule, it has been 
shown [14] that 


Casts 


A / I, /, k(n, t h) t k(ny t }) t k(n, t } 

For H.O and H,S it has been found that &; is essen 
tially zero, whereas k,>k,. The large A found for 
the excited state seems to indicate that ny=1 lf 
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The calculated positions and intensities of the lines are shown below the spectrum 
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state inertial constants presented here may be some- 
What in error due to the inability to assign Y branch 


However, it is felt’ that this 
error is small because of the excellent agrecment 
between the structural parameters of the two mole 
cules \s more bands are analyzed, the ground 
state inertial constants will be continually refined, 
removing any small errors in the present determina 
tion 


transitions uniquely 


The authors express then appreciation to Kk. Bright 
Wilson, Jr. in whose laboratory a part of the analysis 
was carried out 
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p-threose and C'-labeled sodium evanide 





percent yield by means of the lactone 


Synthesis of 0-D-Xylose-1-C“ and B-D-Lyxose-1-C" ' 
H. S. Isbell, H. L. Frush, and N. B. Holt 


Methods are presented for the preparation of a-p-xvlose-1-C"™ and 8-p-lyxose-1-C"™ from 


The xylonie epimer from the evanohydrin svn 
thesis was separated in 34-percent yield by means of the lead salt; the lyxonie epimer, in 49 


Lead p-xylonate-1-C™ was converted quantitatively 


iO D-xvlono-y-lactone-1-C™, and this was reduced by sodium amalgam to a-p-xvlose-1-C™ in 


61 percent \ ield 


vield Thus 


p-Lyxono-1-lactone-1-C" on reduction, gave 8-p-lvxose-1-C" in 64 pereent 
the over-all vields of a-p-xvlose-1-C™ and B-p-lvxose-1-C™ 


based on the C™ 


labeled evanide used in the evanohydrin synthesis, were 21 and 31 percent, respectively. The 


methods reported greatly facilitate the preparation of a-p-xvlose-1-C™ and 8-p-lyxose 


1. Introduction and Discussion 


Polysaccharides D-Xvlose, the most 
i indant of the accompany cellulose in 
wide variety of plant products. For numerous 
vestigations in the fields of bacteriology, plant and 
inimal metabolism, enzymology, and polysaccharide 
chemistry, a supply of 1-C™4labeled xvlose was 
ceeded. Prior to the present study, this substance 

had been obtamed from p-glucose-1-C"™ through the 
niermediate preparation of 1,2 isopropvlidene-p- 
lucose-I-C™, followed by periodate oxidation, and 
eduetion of the resulting 1,2-isopropvlidene-p-rylo- 
trihvdroxyglutaraldehyde [1]. It seemed probable 
that a less complicated svnthesis of D-X\ lose-1-C™ 
vuld be accomplished, and both of the labeled sugars 
-xvlose and p-l\ xose, could be isolated by methods 
similar to those previously emploved at the Bureau 


containing 


pentoses, 


© for the evanohydrin synthesis of other 1-C'-labeled 
> sugars [2]. The method developed is shown in the 
S accompanying outline 

& p- Threose 


: Cyvanohvdrin reaction (NaC'N 
4 * 
CH-OH-COOH (epimeric aldonie acids; p-xvlonie and p 
lvxonic 
: 
a Conversion to lead salts and 
4 separation of epimers 
4 4 
{ p-xvlonate Lead p-lvxonate 
hvdrate not ervstallized 
cation exchange cation exchange 
resil resi 
‘ 4 
i mie acid p-Lyxonie acid 


Lactonizatior Lactonization 


* 
p-Lyxono y-lact one 


»y- lactone 


NaHg, reductior 
. * 


NaHg, reductior 


Separation of salts, 
and ervstallizatior 


Separation of salts, 
and crystallization 
\ vlose-1-C™ 3-p-Lyxose-1-C™ 
1 project on the development of methods f(r the synthesis of radio 
ithohydrates, sponsored by the Atomic Enerey This 
based on the work described in an AEC report 
ke 


t s indicate the literature references at the enc 


Commission 


1 of this paper 
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2. Experimental Procedures 
2.1. 


b-Threose, required as a starting material for the 
cvanohydrin synthesis, was prepared from caleium 
p-xvlonate dihydrate |3] by the Ruff degradation 
14, 5] with certain modifications. Three grams of 
crystalline calcium p-xylonate dihydrate (7.38 mM 
was added to a solution prepared from 0.6 of barium 
acetate monohydrate, 0.15 g of ferric sulfate (approxi 
mately, a hexahvdrate) and 50 ml of water. The 
mixture was heated to 50° C, and treated with 1.8 
ml of 30-percent hydrogen peroxide After the 
solution had darkened, a second 1.8-ml quantity of 
30-percent hydrogen peroxide was added. When the 
solution had again darkened, indicating that reaction 
was complete, it was filtered through decolorizing 
carbon, and concentrated under reduced pressure to 
a thick sirup, which was mixed with 70 mlof methan 
ol, and then with 140 ml of ethanol. The resulting 
precipitate was separated by filtration, washed with 
ethanol, and discarded. The filtrate was concen 
trated under reduced pressure almost to dryness 
and the residual sirup was diluted first with 25 ml of 
5 ml of 2-propanol. The 


»-Threose 


methanol and then with 25 
alcoholic solution was filtered through | g of decoloriz- 
ing carbon, the filter was washed with a 1:1 ethanol- 
2-propanol mixture, and the filtrate was concentrated 
under reduced pressure to remove the alcohol The 
residue was then dissolved in 20 ml of water, and the 
solution was passed through a column containing 
20 ml of mixed cation * and anion * exchange resins; 
the column was washed with 100 ml of water. The 
effluent was concentrated under reduced pressure, anid 
finally was adjusted to a volume of 50 ml. A 5-ml 
aliquot of this solution, when treated with 20 ml of 
N/10 iodine and 30 ml of V/10 sodium hydroxide by 
the method of Kline and Acree [6], reduced iodine 
equivalent to 0.825 mM of aldose. Thus the vield 
of p-threose was 8.25 mM, or 56 percent of the 
theoretical. The material was used without furthes 
purification in the cyanohydrin synthesis 


2.2. Lead p-Xylonate-1-C'* Monohydrate 


A portion of the solution prepared above (20 ml, 
containing 3.3 mM of p-threose) was treated with 








' Amberlite IR 120-H, Rohm & Haas Cx 
‘ Duclite A-4, Chemical Process Co 


Philadelphia, Pa 
Redwood City, Cali 


1.4 mM of sodium bicarbonate, and cooled to 0° C 
It was then mixed with 10 ml of a solution contain- 
ing 2.2 mM of NaC"™N, with an activity of 8.8 me, 
and 2.2 mM of sodium hydroxide.” After 1 day at 
0° C and 4 days at room temperature, the reaction 
mixture was heated at 80° C under a stream of air, 
with the addition of water from time to time, until 
the evolution of ammonia The solution 
was then diluted with water and passed through u 
column of cation exchange resin to convert the salts 
To facilitate removal of carbon 


ceased 


to the free acids 


dioxide, the resin (20 ml) was placed in a tube 25 
mm in diameter, and agitated during the passage of 
the liquid so as to cause the liberated carbon dioxide 


to the surface. The column was washed 
thoroughly with water, and the total effluent, to 
which was added 295 mg of lead carbonate (PbCOs), 
was concentrated under reduced pressure to about 
100 ml, and then filtered. A 500-mg quantity of 
carrier lead p-xvylonate monohydrate was dissolved 
in the filtrate, the solution was concentrated under 
reduced pressure to about 5 ml, and methanol was 
added dropwise to the point of incipient turbidity 
The crystalline lead p-xvlonate-1-C™ monohydrate 
that formed in the course of | day was separated, 
and reerystallized from water by the addition of 
methanol. The two mother liquors were combined 
and concentrated under reduced pressure to remove 
most of the methanol By the successive use of four 
additional 500-mg quantities of nonradioactive car- 
rier, four more crops of lead p-xylonate-1-C"™ mono- 
hydrate were obtained. Finally, the five crops of 
lead p-xvlonate-1-C™ monohydrate were combined 
and recrystallized once from water by the addition 
of methanol. The recrystallized salt weighed 2.557 
g and contained 3,009 ye of activity, corresponding 
to a radiochemical yield of 34.2 percent of the 
NaC'™N used. The residual liquors, containing ap- 
proximately 5,800 we of carbon-14, were used for the 
preparation of p-lyxono-y-lactone-1-C", as described 
in section 2.5 


rise to 


2.3. p-Xylono-y-Lactone-1-C'' 


A water solution of the lead p-xylonate-1-C" de- 


scribed above was passed through a column contain- 
ing 20 ml of cation exchange resin, which was then 
washed with about 1 liter of water. The total 
effluent was concentrated under reduced pressure, 
and the volume was adjusted to 50 ml. Ten-milli- 
liter aliquots of the solution were transferred to 
each of five reduction tubes,® and each tube was 
heated at 60° C under a stream of air until the solu- 
tion had evaporated to dryness. The residue was 
dissolved in a few drops of ethylene glycol mono- 
methyl ether (Methyl Cellosolve), seeded with erys- 
talline p-xylono-y-lactone, and stored at room tem- 
perature over calcium chloride in a desiccator; it 
was moistened each day with a few drops of Methyl 

* A M-percent excess of D-threose was used to insure complete utilization of the 
cyanide 

* Heavy-walled glass tubes, each 200 by 25 mm, with a 24/40 standard-taper 
joint, and an oblique sidearm 50 by 15 mm, attached just below the joint, for 
addition of amalgam after vigorous stirring has begun. A stainless-steel standard- 


taper stopper, with oilless bearings, supports a stainless-steel stirrer, which just 
touches the bottom of the tube 





Cellosolve. After 2 weeks, when the produet 
peared to be entirely crystalline, the material 
used for the preparation of a-p-xylose-1-C" 


2.4. a-p-Xylose-1-C'* 


In order to determine suitable conditions fo: 
duction of the C-labeled lactone, samples of 
radioactive p-xvlono-y-lactone were reduced 
various quantities of sodium amalgam * in the t 
described in footnote 6. The results, given in | 
i, show that the highest vield of reducing sug 
obtained by the use of 4.6 yy of 5-percent son 
amalgam per millimole of lactone. Hence this 
portion was used in the preparation of p-xvlose-| 


m Vv of 


sence of 


Reduction of 1.0 p-rylono y-factor 
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TarpLe | 


sodium amalgam sodium acid o 


nm 


Yield of 
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AV lose 


Sodiurn Naltle, 
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oxalats cent 


Experi 


vent by an 
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* In 2 ml of water 
The five aliquots of p-xvlono-y-lactone-1-C 
described in section 2.3 (each having 600 ye in 1.84 
mM) were reduced by the successive addition t 
each reduction tube of 5.8 g of sodium acid oxalat 
20 ml of ice water, and 8.3 g of 5-percent sodiun 
amalgam in pellet form. The mixtures were stirred 
vigorously in an ice bath until the amalgam was 
spent, at which point the mixtures were combined 
and separated from the mercury. The crystallin 
sodium salts were removed by filtration, washed 
with ice water, and discarded. The filtrate, con 
centrated under reduced pressure to 25 mil, was 
treated at ice temperature with sodium hydroxid 
until it had a faint, but permanent, pink color i 
the presence of phenolphthalein indicator; it was 
then diluted with two volumes of methanol, and th 
resulting precipitate was separated, washed with 
methanol, and discarded after a test showed negli- 
gible radioactivity. A small third crop of salts was 
obtained by reconcentrating the filtrate to 10 ml 
adding 50 ml of methanol, concentrating to a sirup 
and finally mixing the residue with 25 ml of methanol 
The insoluble’ residue, when separated, and 
thoroughly washed with methanol, likewise had no 
appreciable radioactivity. 
he alcoholic filtrate was concentrated 
reduced pressure to remove most of the methanol 
then diluted with water and passed through a column 
containing 20 ml of mixed cation * and anion exchang' 


under 


’ The amalgam was prepared in small pellets by pouring it in molten cor 
through a heated alundum thimble, having a small hole in the bottom 
2ft “shot tower” of mineral ofl. The pellets were stored under oil, « 
before use, were blotted dry, weighed, and rinsed with benzene. 

* Amberlite IR 100-H, Rohm & Haas Co., Philadelphia, Pa. This re 
used in the purification of the labeled sugars in place of the Amberlite I! 
used elsewhere 
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and an additional 5 ml of cation exchange 
at the base of the column. The salt-free and 
ral effluent was concentrated to approximately 
ml under reduced pressure, filtered, and freeze- 
d: the residue was dissolved in 2 ml of methanol 
r 24 hr a crop of crystalline a-p-xylose-1-C™ was 
ined, which, when recrystallized from a few 
yps of water by the addition of methanol, weighed 
mg and had an activity of 1,052 ue. By use of a 
il of 1.1 g of carrier p-xylose in four portions, an 
itional 630 uve of a-p-xvlose-1-C™ was recovered 
m the mother liquor. However, the final mother 
ior contained a substantial quantity of uncrys- 
izable material having an activity of 1,158 ye. 
may have consisted of glycosides formed by 


ition of the sugar with the aleohol used as a 
ent, of oligosaccharides formed by condensation 
wo or more sugar molecules, or of sugar anhy- 
des, together with the polyhydrie aleohol formed 
over-reduction. In order to reclaim p-xylose-1- 
the mother liquor was hvdrolyzed by heating 
Shr with 10 ml of 0.05 N aqueous hydrochloric 
d in a boiling-water bath. The acid was removed 
means of 10 ml of anion exchange resin, and the 
itral effluent was concentrated under reduced 

essure to a sirup. By use of a total of 600 mg of 

‘vlose in three portions as carrier, 154 ye of p- 

vlose-1-C™ was obtained As 1,682 uc of a-p- 

vlose-1-C™ had been separated previously, the total 

or 61 percent of the activity of 
originally used, 


eld was 1.836 ue 
he lead p-xvlonate-1-C" 


2.5. pv-Lyxono-y-Lactone-1-C'* 


solution of the residue from the preparation of 
ead p-xvlonate-1-C' monohydrate (section 2.2 
was decationized by means of a column containing 10 
il of cation exchange resin. One gram of non- 
radioactive p-ly xono-y-lactone was dissolved in the 
efluent The solution was concentrated under 
educed pressure to a sirup and lactonized by heating 
tat 90° C under a jet of air for 8 hr, with the addition 
f 0.2-ml portions of glacial acetic acid at 2-hr 
ntervals. The residue was dissolved in 1 ml of 
ethanol, and 2-propanol was added to the point of 
neipient turbidity. After the addition of seed 
vstals, and storage at room temperature for 24 hr, 
crop of erystalline p-lyxono-y-lactone-1-C™ was 
obtained. When reerystallized once from ethanol 
by the addition of 2-propanol, it weighed 574 mg and 
lad an activity of 2,627 we. By use of an additional 
H00 mg of nonradioactive p-lyxono-y-lactone as 
rier, 1,658 ue of lactone was obtained in three 
portions from the combined mother liquors. The 
tal vield was 4,285 ue, or 48.7 percent, of the 
tivity of the NaC™N used in the original cyano- 
drin synthesis, 


2.6. a-p-Lyxose-1-C'* 


rior to the preparation of a-p-lyxose-1-C™, the 


periments of table 2 were conducted. The pro- 
riions of experiment 3, which gave the highest 


vield of p-lyxose, were used in the preparation of the 
labeled sugar 

In a typical preparation, | mM of p-lyxono-> 
lactone-1-C™, having an activity of 657 we, was 
reduced in the manner described for p-xvlono-y- 
lactone-1-C™, but with the use of 22 ¢ of 4 2-percent 
sodium amalgam,” 8 g of sodium acid oxalate. and 
20 ml of water. Stirring of the preparation was con 
tinued at 0° C for 3 hr, and the resulting 8-p-lyxose 
1-C'* was separated by essentially the same pro 
cedure as that described in section 2.4. A 32-me¢ 
crop of B-p-lyxose-1-C™, having an activity of 140 
ue, Was obtained without carrier, By use of 600 meg 
of carrier p-lyxose, a total of 210 we of the labeled 
sugar was separated in three portions. Hydrolysis 
of the mother liquor, deionization, concentration, 
and use of 0.4 g of carrier in two portions gave 71 ye 
of B-p-lyxose-1-C". The total vield was 421 ue, or 64 
percent of the activ ity of the parent lactone 


Parnes 2 Reduction 


ad m amalgam | 
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